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COMPARATIVE PERFORMANCE ANALYSIS OF VSA AND GWO TECHNIQUES FOR 
LOAD FREQUENCY CONTROL OF SINGLE AREA FUEL CELL MICROGRID 
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Abstract. In this paper, the frequency control of a single area micro grid system is studied. 
Recently, the dynamic performance of micro-grid systems can be enhanced by effective load 
frequency control which is one of the most important aspects. The system having a fuel cell, micro 
turbine, diesel generator, battery energy storage system and flywheel energy storage system is 
modeled. To optimize the dynamic performance of the system, Vortex Search Algorithm (VSA) and 
Grey Wolf Optimizer (GWO) are used. Additionally, a cascade PD-PI controller is also applied to 
the system for enhancing the performance. The results are analyzed with respect to overshoot, 
undershoot and settling time. It can be clearly observed from the results that GWO algorithm provides 
better system performance for selected micro-grid system. 


Keywords: Load frequency control, micro grid, vortex search algorithm, grey wolf optimizer 


СРАВНИТЕЛЬНЫЙ АНАЛИЗ ЭФФЕКТИВНОСТИ МЕТОДОВ VSA И GWO 
ДЛЯ КОНТРОЛЯ ЧАСТОТЫ НАГРУЗКИ ОДНОЗОННОГО ТОПЛИВНОГО 
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В данной статье рассматривается частотное регулирование отдельной микросетевой 
системы. В последнее время динамические характеристики микросетевых систем могут быть 
улучшены за счет эффективного управления частоты нагрузки, что является одним из 
наиболее важных аспектов. 

Смоделирована система, имеющая топливный элемент, микротурбину, дизельный 
генератор, аккумуляторную систему накопления энергии и систему маховичного накопителя 
энергии. Для оптимизации динамических характеристик системы используются алгоритм 
поиска вихрей (VSA) и оптимизатор серого волка (GWO). Кроме того, для повышения 
производительности в системе также применяется каскадный контроллер PD-PI. Результаты 
анализируются относительно времени перерегулирования, недорегулирования и 
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установления. Из результатов ясно видно, что алгоритм GWO обеспечивает лучшую 
производительность системы для выбранной микросетевой системы. 


Ключевые слова: контроль частоты нагрузки, микросетевая система, алгоритм поиска 
вихрей, оптимизатор серого волка 


Introduction: Conventional electrical power systems are mainly consist of large number of 
components. Therefore, transmission of generated power from generation units to the over distance 
last consumer reduces the transmission efficiency [1]. In parallel, these power systems mostly use 
fossil fuels to provide as energy resources. The transfer of energy to long distances causes an increase 
in the fossil fuels used as well. 

In the globalizing world, traditional power systems are gradually being replaced by modern 
power systems [2,3]. Micro-grid systems can be considered as modern power systems. These systems 
usually contain non-fossil source (such as fuel cell, battery systems) and renewable energy sources 
(such as wind turbines, solar panels). 

As with conventional power systems, load frequency control is also a very important 
requirement for micro-grid systems. The main objective of load frequency control can be considered 
as follows; 


1) Keep the system frequency specific band width limits [4], 
2) Deviation of frequency made zero or closest to zero [5], 


In this study, VSA [6] and GWO [7] are applied to the selected micro-grid system to optimize 
load frequency control performances of the system. The frequency deviation obtained by both 
algorithms are compared in terms of overshoot, undershoot and settling times. Moreover, cascade 
PD-PD controller is used in these analyzes. This controller type consists of combination of PD 
controller and PI controller. 

This paper is organized as follows: In section 2, micro-grid model for load frequency control, 
VSA and GWO are briefly mentioned. In section 3, obtained results are illustrated. In section 4, 
Conclusions are given. 


Material and Methods 


Characteristic 
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Fig. 1. Micro grid model for LFC analyzes 


I. Micro grid model for load frequency control: Load frequency control is essential control 
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for micro grid power systems. The micro grid model used in this study is given in Figure 1. This 
model have micro turbine block, diesel generator block, battery energy storage system (BESS) block, 
flywheel energy storage system (FESS) block, micro grid (MG) system block, fuel cell, inverter 
block, connection device block and droop characteristic block. 


In this model, PD controller can be defined as follows: 
PD= Кра + 5Ка (1) 


Where, Краапа ka proportional and derivative gain of PD controller. PI controller can be defined 
as follows: 


ki 


Where, Кр; and ki proportional and integral gain of PI controller. Drop characteristic can be 
defined as follows: 
Drop Characteristic — Z (3) 


Micro turbine can be defined as follows [8]: 


1 


1+STMT (4) 
Where, Tmr is time constant of micro turbine. Diesel can be defined as follows [8]: 
1 
1+5Тр (5) 
Where, Тр is time constant of diesel. Fuel cell can be defined as follows [2]: 
1 
1+5ТЕс (6) 
Where, Trc is time constant of fuel cell. Inverter can be defined as follows [2]: 
1 
1--5Т| (7) 
Where, Ti is time constant of inverter. Connection device can be defined as follows [2]: 
1 
(8) 


1+sTcp 


Where, Tcp is time constant of connection device. MG system can be defined as follows [2]: 
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L (9) 


Ms+D 


Where, M is inertia constant and D is damping coefficient respectively. BESS can be defined 
as follows [8]: 


1 
1+STBESS 


(10) 


Where, Tpess is time constant of battery energy storage system. FESS can be defined as follows 


[2]: 


1 
1+STFESS 


(11) 


Where, Tress is time constant of flywheel energy storage system. In this study, ITSE (integral 
time square error) is used as objective function in GWO and VSA techniques. ITSE can be defined 
as following equation: 


ITSE = f (time x (error)*)dt (12) 
II. Vortex Search Algorithm: VSO is developed by B. Dogan and T. Olmez in 2014 [6]. This 


algorithm has good balance between exploration and exploitation behavior. For this reason, the 
optimal results can be reliably found. The shape of the mixed liquids inspired the development of the 


VSA [6]. 


Generate random solutions 


Tested all candidate solutions 


{ 


Define best iteration solution, reduce the 
radius of circuit, 


i NO 


Made best iteration solution as next circuit 
center 


¥ 


If best iteration solution is better than global 
best, global best= best iteration solution 


Y 


Is termination criteria met? 


YES 


Fig. 2. Flow chart of VSA 
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In this technique, candidate solutions are formed in the first circuit. All candidate solutions are 
controlled for the problem and the best one (iteration best) is determined from these and defined as 
global best solution. Founded iteration best solution is defined as the center of the next circuit. After 
then, radius of the circuit is decreased. New iteration best solution is compared with the previous best 
iteration’s solution. If new iteration best solution has minimum value rather than previous best 
iteration solution, new iteration best solution is defined as global best solution. 

This process continues until the termination criterion is met. Basically flow diagram of VSA is 
illustrated in Figure 2. 


Ш. Grey Wolf Optimizer: GWO is developed by S. Mirjalili et al. in 2014 [7]. This technique 
has very good capability to search of parameter space. Grey wolves are inspired the development of 
the GWO [7]. 


Generate random wolves positions 


Tested all wolves according to the objective 
function 


Define best solution 


If best solution is better than global best, 
global best= best solution 


Update the wolves position 
Is termination criteria met? 
YES 


Fig. 3. Flow chart of GWO 


In this method, positions of grey wolves (candidate solutions) are formed randomly. These 
candidate solutions are search the space for obtained minimum objective function for the problem. 
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After the iteration, best iteration values are determined and then the positions of candidate solution 
are updated. If the best iteration value in the next iteration is better than the previous best iteration 
value, the best iteration value in the next iteration is defined as the global best value. 


This process continues until the termination criterion is met. Basically flow diagram of GWO 
is shown in Figure 3. 


Results and Analysis: In this section, comparative performance results of VSA and GWO are 
discussed. Micro grid system parameters are given in Table 1. In this system, some parameters are 


taken from the reference [2], while some parameters are taken from the reference [8]. 


Table 1: Parameters of micro grid system 


R TMT Tp M D 
0.05 2 2 0.1667 0.015 
Trc Tı Tcp TBEss Tress 
0.26 0.04 0.004 0.1 0.1 


The cascade PD-PI controller parameters with VSA and GWO are given in Table 2: 
Table 2: Controller parameters 


kpi ki kpa ka 
VSA -1.7945 -4.9489 4.9499 2.3019 
GWO 2.8522 4.9500 -4.7296 -1.4668 


Frequency deviations outputs are illustrated in Figure 4: 


x iw 


Frequency Deviation (Hz) 


0 0.5 1 L3 2 2-2 3 
Time (sec.) 
Fig. 4: Comparative frequency deviation graphic 
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From this figure, it is clear that the obtained two frequency deviation outputs with VSA and 
GVVO are close to each other. However, the oscillation amplitude of GWO is slightly larger compared 
to VSA. Obtained numerical overshoot, undershoot and settling time results are given in Table 3. 


Table 3: Comparative performance results 


Overshoot Undershoot Settling time (0.002 band with) 
VSA 3.411e-3 -1.258e-2 0.5054 
GVVO 4.157e-3 -1.236e-2 0.4020 


Itis clearly seen that from this table, although VSA gave minimum overshoot value, GVVO gave 
minimum settling time and undershoot values. 


Conclusions: In this paper, load frequency control of single area having fuel cell micro grid 
system is examined. Comparative performance results performed with VSA and GWO methods. 
Overshoot value obtained by VSA approximately %17.94 lower than GWO. However, GWO gave 
approximately %1.75 lower undershoot and % 20.46 lower settling time value than VSA. It can be 
said that from the obtained results, GWO has the ability to give better performance. 
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THREE-PHASE TVVO-STROKE ELECTROMAGNETIC VIBRATION DEVICE 
MANAGING OUTPUT PARAMETERS 
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Annotation. It is known that the stable operation of existing three-phase low mechanical 
frequency electromagnetic type vibrating devices depends on keeping the specific frequencies of the 
electrical and mechanical systems at the same value or close to each other. In the work process, the 
specific frequency of the mechanical system is variable depending on the mass, while the specific 
frequency of the electrical system remains practically constant. This causes unstable operation of the 
electromechanical system of the vibration device. In order to overcome this deficiency, the article 
considers the issue of expanding the frequency interval of the amplitude-frequency characteristic of 
the electrical system. 

Keywords: low mechanical frequency electromagnetic vibrator, factor of quality, frequency 
characteristic, stable operation, output parameters. 


ÜÇFAZALI İKİTAKTLI ELEKTROMAQNİT VİBRASİYA QURĞUSUNUN 
ÇIXIŞ PARAMETRLƏRİNİN İDARƏ OLUNMASI 


Quliyev H.B., Hüseynov R.A., Yusifov R.A. 
Azərbaycan Texniki Universiteti, huseyngulu @ mail.ru 
Sumqayıt Dövlət Universiteti, ramiz_hu@mail.ru 

“Azərişhıq” ASC, E-mail: ramil77.77@ mail.ru 


Annotasiya. Məlumdur ki, mövcud üçfazalı alçaq mexaniki tezlikli elektromaqnit tipli 
vibrotəsirləndirici qurğular dayanıqlı işi elektrik və mexaniki sistemlərin məxsusi tezliklərinin eyni 
və ya bir-birinə yaxın qiymətdə saxlanmasından asılıdır. İş prosesində mexaniki sistemin məxsusi 
tezliyi kütlədən asılı dəyişən olur, elektrik sistemin məxsusi tezliyi isə praktik olaraq sabit qalır. Bu 
isə vibrasiya qurğusunun elektromexaniki sisteminin dayanıqsız işləməsinə səbəb olur. Bu 
çatışmazlığı aradan qaldırmaq üçün məqalədə elektrik sistemin amplitud-tezlik xarakteristikasının 
tezlik intervalını genişləndirilməsi məsələsinə baxılır. 

Açar sözlər: alçaq mexaniki tezlikli elektromaqnit vibrator, keyfiyyətlilik əmsalı, tezlik 
xarakteristikası, dayanıqlı iş, çıxış parametrləri. 


УПРАВЛЕНИЕ ВЫХОДНЫМИ ПАРАМЕТРАМИ ТРЕХФАЗНОГО 
ДВУХТАКТНОГО ЭЛЕКТРОМАГНИТНОГО ВИБРОУСТРОЙСТВО 


Гулиев Г.Б., Гусейнов Р.А., Юсифов Р.А. 
Азербайджанский технический университет, huseyngulu @ mail.ru 
Sumqayıt Dövlət Universiteti, ramiz_hu@mail.ru 
ОАО «Азеришыг», ramil77.77G mail.ru 


Аннотация. Известно, что устойчивая работа существующих трехфазных 
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низкочастотных вибрационных устройств электромагнитного типа зависит от поддержания 
конкретных частот электрической и механической систем на одном и том же значении или 
близких друг к другу. В процессе работы удельная частота механической системы изменяется 
в зависимости от нагрузку, тогда как удельная частота электрической системы остается 
практически постоянной. Это вызывает нестабильную работу электромеханической системы 
вибрационного устройства. С целью преодоления этого недостатка в статье рассматривается 
вопрос расширения частотного интервала амплитудно-частотной характеристики 
электромеханической системы устройства. 

Ключевые слова: электромагнитный вибратор низкомеханической частоты, 
добротность, частотная характеристика, устойчивая работа, выходные параметры. 


INTRODUCTION 


In many areas of the industry, the application of vibration devices is widely used in the means 
of transportation of finished products, their orientation, mixing, acceleration of chemical reactions 
and other such processes. Since vibration devices are used in all areas of the economy, they have a 
very diverse construction and functional capabilities in a wide power range. The requirement of 
technological processes for vibration devices is that their mechanical frequency is in the required 
interval and stable operation of the system is ensured in this interval.. 

Depending on the values of the mechanical frequency and operating amplitude, their control 
interval, various vibrators are used in technological processes in various fields of industry. The high 
and low values of the mechanical frequency in such devices are very different depending on the place 
of application. There are areas of technological processes in which the use of a low mechanical 
frequency interval (up to 20 Hs) is of great importance [1-4]. Currently, the demand for the creation 
of vibratory devices whose mechanical system operates in the low frequency range is increasing day 
by day. 

It is known that the core of the electromagnets of existing three-phase low mechanical 
frequency electromagnetic vibrators (TTEVD) is assembled from electrotechnical steel sheets. The 
operating principle of the most widespread of such devices is based on the resonance of voltages or 
currents, devices that work with a series or parallel capacitor connected circuit according to the circuit 
of the electromagnet [5,6]. When such devices are operating, the stability of the system depends on 
keeping the specific frequencies of the electrical and mechanical systems at the same or close to each 
other. It is clear that in the work process, the specific frequency of the mechanical system is variable 
depending on the mass, while the specific frequency of the electrical system remains practically 
constant. This causes unstable operation of the electromechanical system of the vibration device. To 
overcome this deficiency, it is more appropriate to expand the frequency interval of the amplitude- 
frequency characteristic of the electrical system. 


1. Theoretical study of output parameters of three-phase two-stroke electromagnetic 
vibration device 


The improvement and renewal of vibration devices used in various fields of industry requires 
the study of the movement of their electromagnetic system. For this purpose, let's study the 
electromagnetic system of the three-phase vibration device, which allows the low mechanical 
frequency oscillation to get moving (figure 1). In the proposed two-stroke three-phase vibration 
device, CA, Св, Cc are the capacities of the capacitors connected in series to the phases, , ra, гв, rc are 
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the active resistances of the phases, La, Lg, Lc are the specific inductances of the phases, Maz, Mca, 


Məc are mutual inductances between phases [1,7]. 
If we write the system of equations of the circuit shown in Figure 1 and show them as a 


system, we get the following: 


Figure 1. The principle diagram of a three-phase two-stroke vibration device 


nz. + Жа m Iade = О. 
JX əl, +126 + Жос =U, (1) 
IjXa + ТАСЫ tZ. m Uc 


Here 
| І | 1 | E 
Z,-7r,Tj я 5 Z=% Tj Dar C. 5 22:11. ar 3 


Хвс = JOM вс; Ха = JOM cy ; X дв = JOM дв. (2) 


It should be noted that the number of windings of the phase windings of the device is equal to 
each other. In addition, the capacities of the capacitors connected in series to the phases are also equal 
in price. So, ra=rB=rc=r; СаА=Св=Сс=С is accepted. 

The specific inductances related to the windings of the phases are obtained from the calculation 
of the magnetic system of the device and are expressed as follows: 


Eie lr 
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L 
м (3) 


-- ————— 
3h + 3uó, 


Here a,,b,,h and а - are the width, thickness, height and distance between the rods of the core, 


respectively. 
And the mutual inductive connections between the coils 


= L; : 
E 2d ” 
+ 
3h +36, 
L 
Мам (4) 
+ — 
3h + 346, 
L L 
Ма = 1,1 = = m 24 


-— 

5 

= 3h 346, 
3h - 3uó 


is expressed in the form. 
According to the given geometric dimensions of the two-stroke three-phase vibration device 


. xo (5) 
3h +36, 
considering that in (4) they are Mas=Mcs=Mca. If we consider inequality (5) in (3) they are 
La=Lp=Lc=L. 

Considering the above in the expressions of total resistances Za, Za, Zc and mutual inductive 
resistances, 


2,-27,-%7с eda ағаш 
ӘС 


(6) 
Xan Xas лан cM cans 
they are. Thus, (1) system equation as we can write 
ire ik ə xəbə kə, 2/7 
LX br 1X) о 0, | (7) 


Бы] =U 
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The low mechanical frequency vibration device consists of a two-stroke three-phase 
electromagnetic core. Such cores belong to the class of ferromagnetic cores, and therefore the main 
property of the material of such a ferromagnetic core is parametrically very similar to the properties 
of existing ferromagnets. At the same time, as a result of the theoretical and practical research 
conducted on the considered two-stroke TTEVD, it was found that if the electrical parameters (Ua, 
Us, Uc and Ia, D», І.) in a two-stroke three-phase system form 120? symmetry with each other, then 
the theoretical-practical operations performed on one of the two tacts can be attributed to both tacts. 
At the same time, since the winding with the same number of windings is used in all three phases, the 
inductances of the phase windings will be the same. However, the magnetic losses due to inter-cycle 
mutual probable periodic current losses are determined by applying the electromagnetic field theory, 
which is directly proportional to the value of the current. Since the electromagnetic losses on both 
cycles vary depending on the length of the air gap in the path of the magnetic flux, the air gap in the 
path of the magnetic flux may differ, even if slightly, in the created two-stroke three-phase vibrating 
device. 

The analysis of the obtained expressions for different phase currents shows that: 

- the time variation of the studied quantities in the steady state is non-sinusoidal in nature, 
which is related to the considerable nonlinearity of the magnetic system; 

- at the beginning of the switching process, the nominal value of the current exceeds the 
nominal value, which is explained by the saturation of the steel of the electromagnets; 

- the non-sinusoidal nature of the oscillations has a great influence on the minimum value of 
the air distance. 

Therefore, it is necessary to ensure the linearity of oscillations in order to facilitate the selection 
of the minimum air distance and to keep its value within reasonable limits. By ensuring the minimum 
value of the air distance, it is possible to obtain a greater electromagnetic pulling force with minimal 
energy losses, which is performed by the correct selection of the parameters of the core, anchor and 
windings of the electromagnet. At the same time, the effect of periodic current loss remains in effect. 
If we solve the cyclic currents generated in the proposed two-stroke TTEVD core using the course of 
the theory of electromagnetic fields, the quantities L and M included in expressions (3) and (4) are 
obtained in complex form: 


L-L,- jL; M - M, - jM,. (8) 


as it is written. Considering these statements in (2): 


jas jo. - п). = (ог, + Да, -Ł), 
ас ас 
(9) 
jx, = jlo(M,j-M,)]=oM, + joM, 
can. If we consider the expression (9) in (5), the following expression is obtained for the work 
performed by the three-phase vibration device in the symmetrical mode: 
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Ir + ob, + Да, Ж +) e" + (oM, + joM ДУ 
ac 


+ (vM, + ем ‚Ге =U,, 


I (OM,  joM e? + 1,(@M, + joM е” + döl 


+(@M ,+ ioM ,)I ,e "^" =U}, 
I (0M , + jM ef +1, (0M, + joM el?" + 


+ (M, + joM ,)I e 7 =U. 


Each equation of system (10) is an independent equation. If the complex form of the phase 
voltages in these equations is expressed by their effective values, then from the first equation of (10) 


3 
p — (11) 
1 А 0 А 0 
r, + oL, + o. nu. + joM ,)e 77 + (oM „+ joM , je?” 
ас 


А 


from the second equation 


2 
гін ау (12) 
r, +@L, + С PR Jen. + ЈМ ,)e 7" + (mM, + дом ,)e 7 
ос, 
from the third equation 
fM _____—_ (13) 
r, + OL, + ІС m J+ loM, + joM , al 4 (ем, + joM, ia 
ос, 
is taken. Here 
jc cp Hc cu 1, v3 
a” =l; a” --—-j—;a-2-— —. 14 
2 и 2 2 , 2 en 


Thus, the denominators of the fractions in the received expressions (11), (12) and (13) are the 
same, and it is enough to study the denominator of one of them to obtain research results for the 
resonance mode of each phase. Let us assume that the frequency of the supply voltage is @ and the 
resonant frequency of the circuit is cə, . Accordingly, voltage resonance occurs in a given series circuit 


and the specific frequency of this circuit 


(15) 
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is vvritten as If vve consider this in the expression of the total resistance of the denominator of 
expression (11), vve get: 


—— 2-2 (16) 
Oy € 
We write from here: 
Z | alL -M о о 
| а а} | 7.9 17 
r, +œ lL, -М,) EA d Mo 


Here, 0, (L, -М 2) is the active resistance transferred from the magnetic circuit to the loop and has 


the same character as the resistance r4. If we make some changes in the last expression and go to 
the module of the expression, we write: 


ЛЕН 0 а | 
dəm a "e 


If we make a transformation in this expression, we write as follows: 


257: (19) 


TA Q €) 


Based on the obtained expression (19), the frequency characteristic of the resistance can be 


established. In this case, the ratio of frequencies -— is accepted as an argument. The term Q, in 
0 


expression (19) was the coefficient related to the current losses is defined as 


"e o,(L, - M,) | (20) 


ra 


The quantity Q, depends on the difference between the active parts of specific inductance and 


mutual inductance and is determined by the expression 


r o,(L, -М,) 


ГА 


0, (21) 


This coefficient is called the quality coefficient of one phase of the three-phase system. 
2. Modeling results of the frequency characteristics of the output parameters of the 
TTEVD 
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The frequency-dependent characteristic of the ratio of the total resistance to the active resistance 


of the phase loop is given in figure 2. The characteristic reflects the dependence of the ratio of total 


resistance (impedance) Z to the active resistance ra of the phase loop on the frequency ratio —. 
M% 


Expression (19) was used to establish these characteristics. Curve 1 shown in Figure 2 is at the value 
of Q, =4 and Q, =3; 2 curve Q, =10, О. =3; 3 curve Q, “100, Q, =3; Curve 4 is constructed 
according to the values of Q, =1000 and Q, =3. 

From the obtained characteristics, it is clear that when the coefficient of quality Q is higher, 
the steepness of the characteristic is obtained more than the curves of Q, =10 and Q, “ 100. Due to 


this, the normal operation of the three-phase vibrator is ensured in the entire frequency range. 
However, the normal operating mode assumed at values of Q, =10, 100 and above is violated [7]. 


If we denote the currents flowing through the circuit of each phase by /,,, corresponding to the 
resonance operating mode and /4 corresponding to the current value of the frequency, we get: 


Lo 1 (22) 


Lio ^ 
(КЕГЕН 


Oy €) 


Based оп the expression (22), the dependence of the modulus of the phase currents on the ratio 


© is established, and these are also shown in figure 3. It can be seen from Figure 3 that when the 
0, 


quality factor changes in the range of 4+100, the parallelism of the resonance curve to the © axis is 
0, 
practically maintained in the specified frequency range of the proposed device. This is due to the 
presence of periodic currents generated in the core of this three-phase vibrating device. It is clear 
from Figure 3 that the rate of frequency-dependent reduction of the amplitude decreases as Q, 
decreases. 
In Figure 4, the {а 


| characteristics are plotted at different values of Q,, with 
AO 


0 
О, =100 remaining constant. It сар be seen from the characteristics obtained here that the steepness 
of the resonance curve decreases due to the increase in the cyclic current losses generated in the core 
of the vibrator. In the figure, curve 1 is constructed at the value of Q, =1, curve 2 at the value of 
О. =3 and curve 4 at the value of Q, =4. 

The same characteristic curves are shovvn for the frequency dependence of the relative 
resistance corresponding to the resonance mode (figure 5). Here, both characteristics are set according 
to the values of Q, =1 and Q, =4. 

Figure 3 shovvs the frequency-dependent characteristics of the phase shift angle at different 
quality factor and constant О, . The changing nature of the curves in this figure repeats the changing 


nature of the curves in figure 4 and figure 5. 
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0.976 0.984 0.992 1.008 1.016 1.024 
Figure 5. Frequency-dependent characteristics of relative resistance 


Thus, the obtained last expressions allovv to determine the dependence of the voltage applied to the 
input of the circuit on the specific inductance and the frequency change. It is known from the course 
of theoretical foundations of electrical engineering that during resonance, the circuit consists of only 
active ( R ) resistance, as the inductive resistance is compensated by the capacitive resistance. This 
frequency obtained in vibration devices is close to the specific frequency of the mechanical system 
and is approximately three times smaller than the frequency of the feeding voltage. On the other 
hand, when 0 > Фо, the frequency change В is positive, and when € < 90, it is negative. When 0) 


=0, В =0 at В = —1 resonance frequency. 
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CONCLUSIONS 


1. The study of the electromagnetic system of the low mechanical frequency TTEVD, 
according to the quality factor corresponding to the resonance mode of its electric circuit, the 
expressions reflecting the complete electrical resistance of the phases, taking into account the current 
losses in the magnetic circuit, and based on these expressions, the characterization of the ratio of the 
current values of the phase currents to its resonance current is assigned. From the obtained 
characteristics, it was found that the steepness of the resonance curve of the system decreases as the 
value of the quality factor decreases. The value of the coefficient of quality Q, =4; When 10 and 100, 


the characteristics obtained in the intended low frequency interval (below 20 Hs) are practically 
parallel to the о/о, line, which makes it reasonable to use P-10 brand steel in the preparation of the 


core. 


2. A mathematical model was obtained with the help of the working parameters and output 
characteristics of the device by designing the equivalent electrical circuit of the low mechanical 
frequency three-phase TTEVD. The theoretically obtained characteristics and report results, which 
correctly and comprehensively reflect the processes taking place in the three-phase vibrator, have 
been confirmed experimentally. 
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MODELS OF ENERGY SYSTEM ELEMENTS IN PHASE 
COORDINATES FOR THE STUDY OF NON-SYMMETRICAL 
STEADY AND TRANSITION REGIMES 


Ibrahimov F.Sh. 
Azerbaijan Technical University 
amfanet @ mail.ru 


Annotation. One of the main issues is the construction of mathematical models of its 
elements and three-phase replacement schemes in order to study the non-symmetric steady-state and 
transition modes of modern energy systems in phase coordinates. Therefore, in the work under 
review, the replacement scheme and mathematical models of the main elements for the calculation 
of the non-symmetric steady-state and transition processes of the energy system are given. 

Keywords: energy system, phase-coordinates, three-phase alternating scheme, non- 
symmetric steady state, transition processes. 


QEYRİ-SİMMETRİK QƏRARLAŞMIŞ VƏ KEÇİD REJIMLORININ 
TƏDQİQİ ÜÇÜN ENERJISISTEMIN ELEMENTLƏRİNİN FAZA 
KOORDİNATLARINDA MODELLƏRİ 


İbrahimov F.Ş. 
Azərbaycan Texniki Universiteti 
amfanet @ mail.ru 


Annotasiya. Müasir enerjisistemlorin qeyri-simmetrik qərarlaşmış və keçid rejimlorini faza 
koordinatlarında tədqiq etmək üçün onun elementlərinin riyazi modellərinin və üçfazalı əvəz 
sxemlərinin qurulması əsas məsələlərdən biridir. Buna görə də baxılan işdə enerjisistemin qeyri- 
simmetrik qərarlaşmış və keçid proseslərinin hesablanması üçün əsas elementlərin əvəz sxemi və 
riyazi modelləri verilmişdir. 

Açar sözlər: enerjisistem, faza kordinatları, üçfazalı əvəz sxemi, qeyri-simmetrik qərarlaşmış 
rejim, keçid prosesləri. 


МОДЕЛИ ЭЛЕМЕНТОВ ЭНЕРГЕТИЧЕСКОЙ СИСТЕМЫ В ФАЗОВЫХ 
КООРДИНАТАХ ДЛЯ ИССЛЕДОВАНИЕ НЕСИММЕТРИЧНЫХ 
УСТАНОВИВШИХСЯ И ПЕРЕХОДНЫХ РЕЖИМОВ 


Ибрагимов Ф.Ш. 
Азербайджанский Технический Университет 
amfanet @ mail.ru 


Аннотация. Одним из основных вопросов современных энергетических систем 
является построение математических моделей ее элементов и трехфазной схемы замещения с 
целью исследования несимметричных установившихся и переходных режимов в фазовых 
координатах. Поэтому в ррассмотренной работе приведены схема замещения и 
математические модели основных элементов для расчета несимметричных установившихся и 
переходных процессов энергосистемы. 
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Ключевые слова: энергетическая система, фазовые координаты, трехфазная схема 
замещения, несимметричный установившийся режим, переходные процессы. 


INTRODUCTION 


Currently, methods and programs for calculating steady-state modes (SSM) and transient 
processes in complex electric power systems (EES) are built on the basis of the direct sequence 
substitution scheme. In order to model transition processes in non-symmetrical short-circuit and 
incomplete phase modes, instead of asymmetry, the option of connecting a shunt or an additional 
resistance was used. In order to determine the additional resistance, it is required to carry out 
additional reports of the phase and line voltage at all branch and node points, taking into account the 
reverse and zero-sequence electrical quantities in the asymmetric mode, and also at the generator 
and load nodes of the circuit in complex networks during the switching process. Therefore, it 
becomes difficult to solve the problem with traditional methods and programs, and it is impossible 
to solve the problem in complex asymmetric modes, when transposition is rejected, when the 
parameters of the line and other elements of the system differ in phase, and also when the three- 
phase transformers feeding distribution substations are connected with a special scheme. 

Taking these into account, it is considered the best method to determine the electrical 
quantities of the circuit in the phase coordinates in steady-state and transient modes [1]. 

The advantage of the a, b, c coordinate system over other systems is that all the mode 
quantities correspond to real values and recalculation of the results is not required to obtain the 
phase quantities. 

The solution of the given problem in phase coordinates is universal to a great extent, as 
various types of short-circuits and phase breaks are easily performed, and the simplicity of the 
algorithm allows it to be used in modern high-performance computers. In this case, the main 
difficulty is to obtain a three-phase replacement scheme for ES elements in non-symmetrical quasi- 
steady and transient modes. Therefore, in the considered work, the replacement scheme and 
mathematical models of the main elements of ES for the calculation of non-symmetric SSM are 
given [1-12]. 


I. MATHEMATICAL MODELS AND SUBSTITUTION SCHEMES IN PHASE 
COORDINATES FOR THE CALCULATION OF NON-SYMMETRIC SSM AND 
TRANSITION MODES OF SYSTEM ELEMENTS 


Mathematical models of the elements of the system and the solution of the problem for the 
description of the replacement schemes in phase coordinates consists of two stages. At the first 
stage, the calculation of the non-symmetric settled mode in the phase coordinates is carried out [7], 
but here it is taken into account that the replacement scheme in the pre-crash mode differs 
noticeably from the scheme of the unilinear mode for the reasons we mentioned above. In the 
second stage, in the adjustable multi-machine ES, the integration of the transition processes in the 
phase coordinates with the reporting program is performed, taking into account the full Park-Gorev 
equations. In the work under review, the equations for the stator loop are formulated using the a, b, 
c coordinate system, and the rotor quantities are formulated using the d, g, 0 system. The periodic 
coefficients in the equations of the synchronous machine are calculated as the angle between the 
stator and rotor axes in each interval of the mathematical solution of the equations. It should be 
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noted that since the periodic coefficients in those equations are expressed by sin y and cosy, their 


calculation does not cause any difficulties. 

In the first stage, the method of calculating the phase coordinates of the three-phase electric 
system was considered. Each three-phase element is illustrated with a corresponding replacement 
diagram. In this case, the operation of power transformers and autotransformers in complex 
transformation ratio for different branching and voltage regulation is taken into account. A matrix 
of nodal equations is constructed for other elements of the system. Based on it, a dialogue complex 
is created that allows to combine the results of the report not only for relay protection and 
automation, but for a more complex issue, i.e. for mode symmetrization, and also with the 
calculation program of electromechanical processes (calculation of short-circuit shunt) in phase 
coordinates. Below is a description of the various three-phase elements of the system in nodal 
equations. Hereafter, all quantities are expressed in relative units (p.u.). 


1.1. A model of a synchronous machine 


If the neutral of the generator is grounded through a conductor as shown in figure 1 and 
considering Í, = /, +1; - 1; - YU, [7] we write: 


р І 
Yi Y» Үз Y, mE U, " 
Y, Y, Үҙ -а ү — 10 U, 
Y. Y Y. aY, =y, Us. 35. 
31 32 = 1 0 3 S Ў | Ü | (1) 
-Y coh -aY 3Y, 0 E, 277 
l "e E, Е, 


here a 217120? = e 77/5 ; У, У,, И — forward, reverse and zero sequence conductors. 


E, 


4 Е, 1. 
z 5 OË A 
n di 6 İ, 


Figure 1. General (a) and generator (b) three-phase replacement schemes of the 
system in case of Y, -Y, 
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The matrix expression (1) includes the total power XS = $, +5, +5, of the synchronous 


generator or induction motor. In the latter case, ES is taken with a minus sign. I, Бала В currents 
are the phase currents at the output of the machine. The in expression (1) the E; emf of phase A, for 
which the calculation is performed, is taken as the main parameter. 

When carrying out calculations in a complex network, one node is taken as a balancing node. 
Then expression (1) for such a node will take the following form: 


Үп Yo Y, E t Ü, 5/0, + HE 

Y Y, Y, 7 С U, E 5,10, алу, Еб’ 

Y, Y, Y, = t Üş 53/0; алу Еб e 
-Y, -Y, -ø 3% +Ү || 0 0 


For a generator unit specified as P, | U |, the total active power is determined taking into 


account Ел. Applying (2), you can obtain the values of apparent and reactive: 


ZS = E Ý (-U, - aU, —aU,)+3Y,| È, | -Enl Un f, 


Q =-1, E, Y, CU, - aU, -a^U,) +3¥,| E, [ - Y Uy Pf 


In this case, the matrix expression (1) can be used as a model of a synchronous machine in 
quasi-deterministic non-symmetrical and short-circuit reporting. 


1.2. Model of power transformers 


To obtain the model of different transformers in phase coordinates, the substitution scheme 
shown in figure 2,a is used without taking into account the magnetization resistance of an ideal 
single-phase two-winding transformer. According to the replacement scheme of the given ideal 
transformer, the replacement scheme of the single-phase two-winding transformer given in figure 
2,b can be obtained. When the nodes k and q are grounded, the replacement circuit turns into a P- 
shaped circuit (Fig. 2,c). 

To adjust the voltage of power transformers (autotransformers) in a complex electrical 
system, it is convenient to calculate the operating modes in its various branches, to recalculate the 
branching resistance of the transformer without changing the resistances of other elements. It should 
be taken into account that recalculation of the resistance in many elements of the system is required 
when branching is connected in the replacement circuit, where the resistances are brought to one 
step of the voltage. In this case, the transformation coefficients of the transformers are recalculated 
to find the true value of the current and voltage. Therefore, it is more appropriate to maintain the 
transformer connection by taking into account the change of branching in the replacement circuits. 

Figure 2,d shows a generalized replacement scheme, taking into account the change of 
branching in the first (а) and second (2) windings of a two-winding single-phase transformer. In 
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MPRA g 


Figure 2. Two-winding transformer and its replacement circuits:a — two 
winding transformer; b - X-shaped substitution scheme, considering the 
transformation factor as (1+t):1; c — P-shaped substitution scheme; d - X- 
shaped substitution scheme, taking into account branching in the first 
(о) and second (/2) circuits 


this case @=1+1t,, B =1+1, here t, and t, shows the change in the number of windings of the 


first and second windings. For example, if the transformation factor is 1.025:1 in p.u., then 
150,025, ір-0. It should be noted that га and tg can be negative and positive. 

In accordance with the generalized replacement scheme of a single-phase two-winding 
transformer, the phase diagrams for different types of transformers and their conductivities were 
obtained taking into account the change of their branching [2]. 

In this case, according to the obtained diagrams, the transformers are shown in the form of 


ГА a matrix of node conductors in phase. 


1.3. Power transmission line model 


The mathematical model of the electric transmission line (ETL) in phase coordinates is 
constructed according to the line model for a single-line circuit. For an example, we will look at 
phase A, which is connected to phase B of a three-phase ETL (Fig. 3, a and b). For the case we are 
looking at, we can write the expression in the form of a matrix below: 
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ав ав ав 
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-Y, їх F ly 
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27 


(3) 


here У, = 27! — it was mutual conductivity between phases a and b; y value of the conductance 


of the shunt connected to ground; 0,0, və Tor — are the node voltages and currents, 


respectively. 


(3) for the three-line model, taking into account the transverse interconnections and 
conductances to ground (shunt) of the ETL, can be vvritten according to Fig. 3, c and d: 
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вс |’ (5) 


In this case, comparing the one-line diagram of ETL with the three-line diagram, it can be 


seen that each element in the three-line consists of a 3x3 matrix. 
2% (Or Y 


aec 


) resistance matrix (5) can be determined according to the basic parameters in a 


suitable straight way or with the help of a transformation matrix including pre-calculated symmetric 
complexes (Zo, Zi, Z2). If the matrix of symmetric complexities is defined as a diagonal matrix 


Zo12, then we can write: 


| ТАРТИ 


е | п 
|= 3| T 1 Zo IT | 5 


where | Т li the transformation matrix is determined: 


1 1 d 
= а’ a 
la a 


(6) 
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Here, in the expression (6), it is taken into account that the straight and reverse sequence 
resistances (om/km, f=50 Hs) of the three-phase single-cycle ETL are the same. Therefore, the 
matrix Zaee determined according to expression (6) is symmetric. 

Taking into account the above, expression (5) can be written as follows: 


Zaa ав Zac 
Z aec = Zae із Zec d (7) 
Lac Zec Zec 


Figure 3. Scheme of three-phase ETL 
a- single line diagram of three-phase ETL; b- replacement scheme of ETL, taking into account the 
conductivity to ground; c- replacement scheme of the ETL, taking into account its transverse 
connection; d- replacement scheme of the ETL, taking into account the transverse connection and 
conductances (shunts). 


In other words, the number of (N) elements forming the matrix (7) is 6. In order to determine 
the elements of the resistance (conductance) matrix for three-phase circuits, including single-line 
circuits, the following expression is used: 
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N= 1 (n+n’) 
2 
here, n — is the total number of wires. For example, for a single-phase line N=/, and for a three- 


phase line N=3. If the fourth wire of the three-phase ETL is also taken into account (figure 4), then 
N=10. Therefore, the symmetric matrix Z,,. will have the following form: 


Z ГАЖ Za s 
2 ав Zee Loc ЖА 
ed Za Zec 7. а (8) 
Z Ze Zo Zə 
Figure 4. Four-vvire scheme of three-phase ETL 
Considering expression (8), the knot matrix expression takes the form belovv: 
. . 1. А . 
IL Y us + 2 Үлш + I- ^ d aec3 U, 
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1.4. Load model 


In the reports, the unsymmetrical load can be given according to the value of the full power 
on the phases. In this case, if the adjustment ranges of the transformers are defined clearly enough, 
then they can be given in the form of $ н =P, + jQ, = const in their non-symmetric reports. It 
should be noted that the value of voltage is not the same in case of asymmetrical distribution of 
loads according to phases. Therefore, in the calculation of non-symmetrical electrical SSM, the load 
Y,, =const cannot be considered as a shunt with constant conductance. The value of Ү,, is 


нш 


determined separately for each phase. In this setting, mode reporting can only be done with iteration 
methods. 


П. MODELING RESULTS FOR SYMMETRIC MODES IN PHASE COORDINATES 


Table 1 
Parameters of the generator 
Type of machine salient pole 
The voltage 6kV 
Power 1290 kVA 
Speed of rotation 300 cycle/min. 
Number of poles 20 
The voltage of the exciter 80 V 
Resistances X, | 0,957 p.u. 
X; | 0,345 р.а. 
X; | 0,193 р.а. 
X, | 0,196 p.u. 
X; | 0,282 p.u. 
Хо | 0,0573 p.u. 
Time constants Т, | 1,03 $ 
Ту | 0,052 s 
T, | 0,0816 s 
Parameters of the exciter 
Xy 0.17 p.u. 
X, | 0.15 p.u. 
R, | 0,2 Om, 75° 
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Calculation experiments were carried out for the synchronous machine test sample scheme 
[15] working in the system for the purpose of visual representation of the developed methodology 
for the description of the energy system elements in phase coordinates and the software's 
functionality. Initially, the symmetric mode was considered. All the necessary data of the generator 
for the studied system are given in table 1. 


The report results in the phase coordinates for the three-phase short-circuit at the generator 
outputs are given in table 2 in the mode with power flows, and in table 3 in the case of power flows. 
The obtained table results and the experimental curves of the transition process (Fig. 5) completely 
coincide, which confirms the adequacy of the models of system elements presented in phase 
coordinates. 


Ig р.и. 


4 


Figure 5. Variation curves of the induced current (a) and the 
current of phase A (b) during a three-phase short circuit 
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Table 2 


Results of calculation of SSM in phase coordinates (simmetric) 

Nod | Node iU | Angle P, Qe РІ 01 
eno | type (ou) (deg) (MW) (MVAR) (MW) (MVAR) 
1 3 1,056 -1,527 0,000 0,000 0,000 0,000 
2 3 1,056 -121,656 | 0,000 0,000 0,000 0,000 
3 3 1,056 118,479 0,000 0,000 0,000 0,000 
4 3 1,054 -1,540 0,000 0,000 0,000 0,000 
5 3 1,054 | -121,575 0,000 0,000 0,000 0,000 
6 3 1,054 116,490 0,000 0,000 0,000 0,000 
7 3 1,055 -2,431 0,000 0,000 0,000 0,000 
8 3 1,055 -122,460 | 0,000 0,000 0,000 0,000 
9 3 1,055 117,608 0,000 0,000 0,000 0,000 
10 3 1,048 87,473 0,000 0,000 8,333 5,000 
11 3 1,048 -32,453 0,000 0,000 8,333 5,000 
12 3 1,048 -152,507 0,000 0,000 8,333 5,000 
13 3 1,043 87,153 0,000 0,000 10,000 6,333 
14 3 1,043 -32,757 0,000 0,000 10,000 6,333 
15 3 1,043 -152,809 | 0,000 0,000 10,000 6,333 
16 3 1,045 87,183 0,000 0,000 11,670 7,333 
17 3 1,045 -32,741 0,000 0,000 11,670 7,333 
18 3 1,045 -152,736 | 0,000 0,000 11,670 7,333 
19 3 1,043 87,099 0,000 0,000 10,000 6,333 
20 3 1,043 -32,819 0,000 0,000 10,000 6,333 
21 2 1,043 -152,861 0,000 0,000 10,000 6,333 
22 2 1,043 87,086 0,000 0,000 13,333 8,333 
23 2 1,043 -32,830 0,000 0,000 13,333 8,333 
24 3 1,043 -152,872 | 0,000 0,000 13,333 8,333 
25 1 1,070 0,000 32,373 16,230 0,000 0,000 
26 1 1,070 | -120,072 | 32,306 15,976 0,000 0,000 
22 1 1,070 120,072 | 32,497 16,081 0,000 0,000 
28 5 0,000 0,000 0,000 0,000 0,000 0,000 
29 4 1,070 0,185 19,985 10,151 0,000 0,000 
30 4 1,070 | -119,839 | 20,105 10,292 0,000 0,000 
31 4 1,070 120,209 | 19,912 10,329 0,000 0,000 
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32 5 0,000 0,000 0,000 0,000 0,000 0,000 
39 3 1,057 -2,385 0,000 0,000 6,000 2,333 
34 3 1,057 -122,414 | 0,000 0,000 6,000 2,333 
35 3 1,057 117,649 0,000 0,000 6,000 2, 333 
36 4 1,070 -1,522 10,047 8,272 0,000 0,000 
37 4 1,070 | -121,546 | 10,075 8,292 0,000 0,000 
38 4 1,070 118,502 9,880 8,428 0,000 0,000 
39 5 0,000 0,000 0,000 0,000 0,000 0,000 
40 3 0,521 -93,090 0,000 0,000 0,000 0,000 
41 3 0,855 83,150 0,000 0,000 0,000 0,000 
42 3 0,836 -5,011 0,000 0,000 0,000 0,000 
43 3 0,297 -76,996 0,000 0,000 0,000 0,000 
44 3 0,004 | -175,166 | 0,000 0,000 0,000 0,000 
45 3 0,602 67,101 0,000 0,000 0,000 0,000 
46 3 0,602 -32,817 0,000 0,000 0,000 0,000 
47 3 0,602 | -152,660 | 0,000 0,000 0,000 0,000 
48 3 1,043 -122,888 0,000 0,000 0,000 0,000 
49 3 1,043 117,158 0,000 0,000 0,000 0,000 
50 3 1,043 -2,545 0,000 0,000 0,000 0,000 
Table 3 
Reporting results of SSM in phase coordinates (symmetrical mode, power flows) 
Start Endof | P (MW) Q Start of End of P Q 
of node (MVAR) node node (MW) | (MVAR) 
node 
1 23 -32,373 | -15,122 25 1 32,373 16,230 
2 26 -32,306 | -14,879 26 2 32,306 15,976 
3 27 -32,497 | -14,971 27 3 32,497 16,081 
4 29 -19,985 -9,405 29 4 19,985 10,151 
5 30 -20,105 -9,534 30 5 20,105 10,292 
6 31 -19,912 -9,582 31 6 19,912 10,329 
10 16 12,079 3,889 16 10 -12,063 -4,690 
11 17 12,106 3,965 17 11 -12,090 -4,768 
12 18 11,943 4,030 18 12 -11,926 -4,833 
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10 22 12,317 5,404 22 10 -12,292 -6,470 
11 23 12,230 5,382 23 11 -12,207 -6,452 
12 24 12,103 5,536 24 12 -12,077 -6,610 
13 19 6,511 1,803 19 13 -6,509 -2,090 
14 20 6,437 1,662 20 14 -6,436 -1,950 
15 21 6,466 1,749 21 15 -6,465 -2,037 
13 22 4,002 0,986 22 13 -4,001 -1,568 
14 23 3,912 0,880 23 14 -3,911 -1,463 
15 24 3,974 0,972 24 15 -3,973 -1,555 
16 19 4,617 3,748 19 16 -4,614 -4,470 
17 20 4,489 3,649 20 17 -4,487 -4,373 
18 21 4,444 3,785 21 18 -4,440 -4,508 
19 22 0,996 0,186 22 19 -0,996 -0,626 
20 23 0,924 0,139 23 20 -0,924 -0,579 
21 24 0,988 0,203 24 21 -0,988 -0,643 
7 33 -3,966 -6,288 33 7 3,969 5,701 
8 34 -4,024 -6,332 34 8 4,028 5,745 
9 35 -3,837 -6,442 35 9 3,841 5,856 
33 36 -10,047 -8,021 36 33 10,047 8,272 
34 37 -10,074 -8,038 37 34 10,075 8,292 
35 38 -9,881 -8,178 38 35 9,880 8,428 
CONCLUSIONS 


1. For the investigation of non-symmetrical steady-state and transition modes of complex 
energy systems, phase coordinate models of system elements (synchronous generator, transformer, 
power transmission line, load) were developed and the algorithm of using these models in the 
software was given. It has been shown that when reporting non-symmetrical quasi-decided modes 
in phase coordinates, it is not appropriate to accept the load, unlike the unilinear schemes. When the 
loads on the phases are unevenly distributed, the voltages on them are not the same. Therefore, 
during such reports, it is suggested to consider loads in the form of constant conductance of each 
phase. 

2. The proposed models of system elements can be used for the study of any type of 
longitudinal-transverse damage, complex non-symmetrical quasi-steady regimes, short-circuit 
transition processes (without preliminary calculation of short-circuit shunts) in phase coordinates. 
The dialogue software complex developed on the basis of the proposed models can be applied 
during the management of modes in operative-dispatcher management in the case of any type of 
automation. 
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3. Reports were made for the three-phase short-circuit mode in the phase coordinates of the 


synchronous generator scheme operating in the electric network for the cases without and with load 
currents. The coincidence of the change curves of experimentally obtained currents and the reported 
results shows the adequacy of the proposed models, the methodology and algorithm developed on 


its basis. 
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Abstract: Electricity as a commodity has unique characteristics. So, since it is not possible to 
store it in large volumes, the buying and selling processes cause more difficulties compared to other 
commodities. Delivery of electricity from production to consumption in a very short time (tenth of a 
second), transmission and distribution networks having a complex infrastructure is one of the main 
reasons for the inevitable technical losses of electricity. Available the problem does not end there. If 
electricity is not metered during its short journey from production to consumer, it causes additional 
losses. Due to the improper use of electricity and the ineffective management of transmission and 
distribution networks, not only technical losses occur at the stages from production to consumption, 
but also non-technical, i.e., commercial losses. Despite the fact that a lot of work has been done to 
seriously reduce these losses in recent years , this is still the case remains an important problem. 
Thus, in 2021, 575.0 million kV of unbalance, i.e. losses exceeding the normative technical loss, 
occurred in the regional networks of "Azerishiq" OJSC. The conducted studies show that reducing 
these losses to the level of normative technical losses, i.e. up to 8.8%, is economically more 
efficient than putting new production capacities into use. 

At a time when energy security is the guarantor of the independence of countries, in 
connection with the transition to a liberal market economy, bringing the losses to the minimum 
level during the transmission and distribution of electricity is one of the most important tasks facing 
the states. 


Keywords: distribution networks, energy balance, actual losses, technical losses, commercial 
losses, load losses, conditional-stable losses, climatic losses. 


ELEKTRİK ŞƏBƏKƏLƏRİNDƏ ELEKTRİK ENERJISININ BALANSI VƏ 
FAKTİKİ İTKİLƏRİN TƏHLİLİ 


ÇİNGİZ ŞƏFİ OĞLU İBADOV 
iqtisad üzrə fəlsəfə doktoru 
“Azərişıq”Açıq Səhmdar Cəmiyyəti 

e-mail: chingiz.ibadov @ mail.ru 


Annotasiya: Elektrik enerjisi bir əmtəə olaraq unikal xüsusiyyətlərə malikdir. Belə ki, onun 
böyük həcmlərdə saxlanılması, yəni anbarlaşdırılması mümkün olmadığından alış-satış prosesləri 
digər əmtəələrə nisbətən daha çox çətinliklər yaradır. Çox qısa bir vaxt ərzində (saniyənin onda biri) 
elektrik enerjisinin istehsaldan istehlaka qədər çatdırılması, ötürücü və paylayıcı şəbəkələrin 
mürəkkəb infrastruktura malik olması elektrik enerjisinin texniki itkilərə məruz qalmasını 
qaçınılmaz edən əsas səbəblərdən biridir. Mövcud problem bununla da bitmir. Elektrik enerjisi 
istehsaldan istehlakçıya qədər qısa zaman kəsiyində keçdiyi yolda ölçülməzsə, bu da əlavə itkilərə 


EEEC Volume 13 Nel 2023 37 


səbəb olur. Elektrik enerfisindən düzgün istifadə edilməməsi, ötürücü və paylayıcı şəbəkələrin 
səmərəli idarə edilməməsi səbəbindən istehsaldan istehlaka qədər olan mərhələlərdə nəinki texniki 
itkilər, hətta əlavə olaraq qeyri-texniki, yəni kommersiya itkiləri də baş verir. Bu itkilərin son illərdə 
ciddi olaraq azaldılması istiqamətində bir çox işlərin görülməsinə baxmayaraq, bu hələ də mühüm 
bir problem olaraq qalır. Belə ki, “Azərişıq” ASC-nin region şəbəkələrində 2021-ci ildə 575,0 mln 
kVs qeyri-balans, yəni normativ texniki itkidən artıq itkilər baş vermişdir. Aparılan araşdırmalar 
göstərir ki, bu itkilərin normativ texniki itkilər səviyyəsinə, yəni 8,8%-э qədər endirilməsi yeni 
istehsal güclərinin istifadəyə verilməsindən iqtisadi cəhətdən daha səmərəlidir. 

Enerji təhlükəsizliyinin ölkələrin müstəqilliyinin təminatçısı olduğu bir dövrdə, liberal bazar 
iqtisadiyyatına keçidlə əlaqədar, elektrik enerjisinin ötürülməsi və paylanması zamanı itkilərin 
minimum səviyyəsinə çatdırılması dövlətlər qarşısında duran ən mühüm vəzifələrdən biridir. 


Açar sözlər: paylayıcı elektrik şəbəkələri, enerji balansı, faktiki itkilər, texniki itkilər, 
kommersiya itkiləri, yük itkiləri, şərti-stabil itkilər, klimatik itkilər 


As a result of the reforms carried out under the leadership and initiative of the President of the 
Republic of Azerbaijan Ilham Aliyev, in accordance with the strategy of the founder of our 
independent state, National Leader Heydar Aliyev, aimed at the dynamic development of the 
country's economy, the industrial and socio-economic potential of our republic has been further 
developed and diversified, and high achievements have been achieved in all areas of the economy. 
the economic security of our sovereign state has been ensured. 


The recent global economic crisis has had an impact on the electricity sector as compared to 
other sectors. It is for this reason that the constant attention of the President of the country directed 
to the electric power sector as a whole, the implemented reforms, especially in the territories 
liberated from occupation, which have great economic potential and productive conditions. The 
successful implementation of the "Smart City", "Smart Village" and "Green Energy" projects based 
on the application of advanced innovations, modern and digital technologies will allow Karabakh 
economic region to become one of the most developed regions of our country in a short time. As a 
result of the rapid socio-economic development in the country, new power centers are being created 
to meet the electricity demand of newly created industrial enterprises, cultural and social objects, 
agriculture and other infrastructures, distribution networks are being rebuilt based on new 
technologies and in accordance with international standards, circular electricity the number of 
networks is increased, immediate and necessary measures are taken in the direction of efficient 
delivery of electricity to each consumer. Legislation in the field of electric energy is constantly 
being improved, structural and personnel changes are being made. 


The government's attention to transmission and distribution networks is constantly 
increasing. Reconstruction works are proceeding rapidly. However, despite this, the total volume of 
losses in our country (commercial and technical losses) is still higher than the level of the 
corresponding indicator of developed countries. In European countries, this rate is 6%. Therefore, 
the amount of electricity consumed must be monitored, measured and calculated. One of the most 
important methods used for this process is the balance method. The electric energy balance is the 
basis for planning and optimizing power plant operating modes, fuel reserves, equipment 
maintenance schedules, electricity transmission, distribution and sales service tariffs. 
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In the electric network means the equality of the electric energy supplied to the electric 
network (wholesale and retail market subject) and the electric energy sold to consumers, taking into 
account the actual losses for the transmission and distribution of electric energy and the production 
needs of the network organization. 


The amount of actual losses is defined as the difference between the amount of electricity 
supplied to the grid from other networks or electricity producers and the amount of electricity 
consumed by subscribers connected to the grid, taking into account the electricity transmitted to 
other grid companies. 


The actual losses for the transmission and distribution of electricity through networks are the 
most important indicator of the economic efficiency of the operation of networks, an indicator of 
the state of the entire energy system, and a clear indicator of the interrelationship of the efficiency 
of energy sales activities. 


Liberalization of the electric power system and utility sector are carried out in our country, 
the creation of a large number of structural units, companies, distributors, regional and wholesale 
companies, depending on the technology and types of work, makes the task of forming the 
structural components of the balance of actual losses even more urgent. Thus, it significantly 
increases the requirements for the reliability of the balance statements, it requires the preparation 
and approval of normative-legal documents regulating the procedure for determining the balance in 
electric networks, as well as the procedure for calculating its structural components. The importance 
of solving these problems and the complexity of the work can be explained as follows. When the 
distribution network receives electricity from wholesale, it is obliged to compensate (pay) taking 
into account the actual losses of electricity in its network. In other words, the distribution network is 
obliged to pay for the energy it buys but cannot sell, that is, it actually loses it. Therefore, the value 
of normative technological losses of electric energy should be consider and taken into account 
while calculating the tariff for electricity transmission, distribution and sale services. If the actual 
losses in the distribution and transmission network companies are greater than the normative 
technological losses, economic efficiency cannot be achieved in these networks. Thus, distribution 
and transmission networks pay for all actual losses, receive compensation in the amount of the 
value of normative technological losses of electricity (taken into account in the tariff), endure 
financial losses due to "excessive" losses of electricity and operate at a loss . The "excess" of 
electricity losses is due to errors in the metering system, incomplete accounting of consumed 
electricity, cases of theft, low ability of consumers to pay, ineffective control mechanisms to reduce 
electricity losses in networks and many other reasons. [1] 

That is why the main task of the transmission and distribution networks is to accurately 
determine the electrical energy balance and the structural components of the balance, including the 
information about the consumed electrical energy, i.e. to increase the validity and reliability of the 
calculation of the actual losses. 


In order to determine the role of balance sheets and the amount of actual losses in the 
energy systems of countries, let's consider the indicators of our republic, the Russian Federation and 
some former Soviet Republics, as well as several foreign countries. Table 1 shows the dynamics of 
electricity losses in the power grids of the Russian Federation [11]. 
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Table 1 
The dynamics of electricity losses in the power grids of the Russian Federation 
Unit of 
Balance 
A measur | 1994 | 1995 | 2000 | 2004 | 2005 | 2008 | 2010 | 2015 | 2020 
indicators 
ement 
Produced billion 1040. | 1037, | 1049 | 1089. 
1. əl ” 
electricity kW —... 38 7 .9 67 
Specific 
нат аса астаса ағы 
n of power kW 1 
plants 
Purchased billion 
electricity kW 222 | 184 | 102 12.2 10.1 3.11 2.1 8.2 | 1,378 
(import) 
Electricity DM 
: billion 
supplied kW 41.7 38 22.8 19.8 22.5 | 20.74 | 19.5 | 182 | 12,11 
(export) 
Released to | billion 781. 955.3 972. | 1009. 
desean kW 794.7 8 803.5 | 864.9 | 875.7 5 956.1 89 35 
Consumed | billion 698. 846.1 859. 
707. 701.9 | 752.2 | 763.1 1.1 l 
electricity kW we 3 “. ui 0 р 33 də 
billion 109.2 | 104.9 | 113. 
79 83.5 | 101.6 | 112.6 | 112.6 103.5 
Electricity kW 4 7 5 
ii % | 994 S 12.64 | 13.02 | 12.86 | 1143 | 10.98 P 10.25 


It can be seen from Table 1 that, while the electricity supplied to the network increased by 
7.19% for the period 1994-2004, the losses increased by 37.64% in 10 years. Since 2005, a 
downward trend has been observed. 


Now let's analyze the information about the electricity consumption and losses of several 
former Soviet Republics, which I obtained from the official website of the "International Energy 
Agency" (Table 2). 


Table 2 
Data on consumed electricity and losses in CIS countries 
2000 2005 2010 2015 2020 
Consu Cons Consu Consu 
med umed Gülən med med 
DES Loss anei Loss ed Loss не Loss “ürey: Loss, 
, 70 , о | energy, | ,% , 70 % 
n y, M M Watts a Ын 

Watts Watts Watts Watts 

: 18.0 | 175.3 | 13,3 11.4 | 152,34 133,49 
Ukraine | 159,955 8 5 4 178,313 9 8 Tel 7 10.4 
Belarus | 24,537 js gp l d 32.81 P 32,056 | 12.1 | 35,936 | 11.4 
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Kazakhs 48,621 13.9 | 64,16 | 10,2 78.09 8.03 | 86.696 14.0 | 103,59 18.5 
tan 6 4 1 7 
16. 15. 10. 
Georgia 77 i 7,143 i 9,985 с. 10,605 | 6.6 | 11,609 T3 
5 44,323 | 9.08 d 8.84 | 49,266 | 8.82 | 54,552 | 8.8 | 61,046 8.8 
251” ы e ^^ оог 2” 
tan 9 8 9 
ЕНЕГЕ 16.92 ЕЛ 6 Er, ИС 
n 4 7 8 0 
Table 3 
Information about consumed electricity and losses in foreign countries 
2010 2015 2020 
С 
- Consum Consum 
TN Loss, | Loss ed Loss, | Loss, ed Loss, | Loss, 
ти M Vts | ,% | energy, | М Vts % energy, | M Vts % 
? 
Watts M Watts M Watts 
Belgium 90251 | 4283 | 4.75 | 86917 3 815 4.39 83 117 3 444 4.14 
Austria 65 422 | 3350 | 5.12 | 67050 3 466 5.17 66 911 3191 4.77 
Czechia 62 559 | 4466 | 7,14 | 62 740 4 066 6.48 63 271 4117 6.51 
Denmark | 35737 | 2624 | 73341 33832 1 765 2:22 34 194 1573 4.60 
Finland 87 720 | 2761 | 3.15 | 82490 2434 2.95 81 869 3021 3.69 
France . 35 414 | 6.99 | 484187 | 36140 | 7.46 | 459578 | 35893 | 7.81 
570 
Germany 840 23 973 | 4.20 553955 | 25605 | 4.62 525180 | 26939 | 5.13 
Greece 59058 | 3782 | 640 | 57340 4 895 8.54 52 315 3 256 6.22 
Hungary 39807 | 3800 | 9.55 | 41836 3 695 8.83 44 591 3 139 7.04 
Ireland 27 622 | 2098 | 7.60 | 27768 2 084 7.51 30 832 2313 7.50 
Italy . 20 570 | 6.22 | 316 897 | 19717 | 6.22 | 301755 | 17760 | 5.89 
Estonia 8 477 1047 | 123 8 136 697 8.57 8 588 404 4.70 
Latvia 6 942 725 10.4 6914 450 6.51 7 009 381 5.44 
Lithuania | 10 316 989 9.59 | 10,960 793 7.24 11,975 951 7.94 
. 117 
Netherlan 600 5 633 | 4.79 | 114651 | 5264 4.59 | 116491 | 5059 4.34 
ds 
130 
Norway 100 9554 | 7.34 | 127647 | 7469 5.85 | 133180 | 8909 6.69 
141 
Poland 275 11 851 | 8.39 | 149461 | 10534 | 7.05 | 158 080 9995 6.32 
Portugal 54 893 | 4280 | 7.80 | 51 746 4 894 9.46 51 551 4501 8.73 
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Slovakia | 25936 | 856 МЫ 26737 | 1363 ER 23341 | 1571 ERE 
Slovenia | 13045 | 982 | 7.53| 13758 | 864 | 628 | 11360 | 849 | 747 
Spain 2 27 400 | 9.85 | 265 098 | 26 509 | 10.0 | 262 409 | 25631 | 9.77 
147 
Svreden io, | 10561 | 7.17 | 136263 | 6988 | 5.13 | 135287 | 10434 | 7/71 
E 64157 | 4372 | 681 | 62787 | 4540 | 723 | 59894 | 4190 | 7.00 
364 
England 395 | 26884 | 738 | 371954 | 28719 | 7.72 | 313760 | 26323 | 8.39 
202 
Turkey >> | 30222 | 149 | 253841 | 36528 | 143 | 292070 | 29126 | 997 
, 236 
Australia | %7, | 16381 | 6.93 | 234756 | 14296 | 6.09 | 250020 | 8933 | 3.57 
New 
43446 | 3101 | 7,14 | 42895 | 2902 | 636 | 43768 | 2755 | 629 
Zealand 
559 
Canada jg, | 48659 | 8.69 | 578 723 | 33263 | 575 | 565814 | 33650 | 5.95 
: 1155 | 2 4 150 255 4057 197 
America 366 999 6.28 687 322 6.15 937 969 4.88 
Mexico = 44252 | 16.8 | 297107 | 40640 | 13.6 | 334805 | 31633 | 945 
416 
Korea тә | 18034 | 3.78 | 524141 | 17979 | 343 | 557463 | 18433 | 3.31 
1102 1025 
Japan oog 47917 |434 |) | 41162 | 401 | 993136 | 44194 | 445 


The loss percentages shown in this table are not the actual losses, but rather the indicators of 
the normative technical losses (technological consumption) used in all former Soviet countries 
(approved by the Cabinet of Ministers of the Republics). The actual losses are not calculated, the 
balance of the energy system is drawn up formally. 


In table 3 below, let's consider the tables about the consumption and losses of electricity 
prepared by the International Energy Agency based on the data provided by foreign countries [12]. 


In Austria, the Czech Republic, Finland, France, Norway, Portugal and Sweden, electricity 
losses are determined by the balance method as the difference between electricity delivered to the 
grid and consumption, i.e. actual losses. 


In the table 4-6 below, let's consider the comparative tables of "Azerishiq" JSC for the years 
2016-2021 on the purchased, sold and actual losses [14]. 


Table 4 
Comparison of actual losses of "Azerishiq" OJSC Baku RETSI for the years 2016 - 2021 
Year | Purchased | Total Sold 1000 Commercial 
EA НГ Actual losses technical 
5 Electricity | Electricity Р losses 
requirements 
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kW kW kW % kW % kW % 
9 220 282 8 771 638 | 448 644 673 205 -224 561 | -2.44 
-. 598 200 398 TRA 889 -— 491 % 
8 430 181 8 024 691 | 405 490 631 658 -226 167 | - 2.68 
— 864 120 744 Siue 089 TARI 344 % 
8 978 613 8 496 744 | 481 869 708 007 -226 138 | - 2.52 
.. 450 331 119 — 231 — 112 % 
8 953 873 8 338 920 | 614952 714 485 -99 532 | - 1.11 
ə 311 207 715 -— 575 -— 860 ZA 
8 739 055 8 163 777 1 575277 703 059 -127 781 | - 1.46 
xn 170 311 860 PNE 507 -— 647 % 
8 963 210 8 221 704 1 741 505 737 086 4419 0.05 
— 418 720 698 а 374 — 324 % 
Table 5 
Comparison of actual losses in the regions of "Azerishiq" OJSC for the years 2016 - 2021 
Purchased | Total Sold Orme Commercial 
Year NE .. Actual losses technical 
Electricity | Electricity Е losses 
S requirements 
kW kW kW % kW % kW % 
10542 648 | 8911 156 1631491 | 15.4 | 1058 462 | 10.0 5.4 
oer 155 677 478 8 211 4 — 4 
10 104 858 | 8427 708 1677 150 | 16.6 | 1 025 037 | 10,1 6.4 
- 601 065 535 0 138 4 —. 5 
10 056 771 | 8 341616 1715 154 | 17.0 | 1131268 | 11.2 5.8 
55 364 874 489 2 031 5 а 1 
9 831 736 7 978 805 1852930 | 18.8 | 1 117088 | 11.3 7.4 
р 544 867 677 2 355 6 — 8 
2017 9 848 611 7 600 048 2 248 563 | 22.8 | 1128300 | 11.4 | 1120262 | 11. 
369 088 280 3 408 6 872 37 
2016 10 293 253 | 7824 095 2 469 157 | 23.9 | 1196923 | 11.6 | 1272234 | 12. 
174 646 528 9 069 3 458 36 
Table 6 
Comparison of actual losses in the networks of "Azerishiq" OJSC for the years 2016 - 2021 
Purchased | Total Sold 1. Commercial 
Year n E Actual losses technical 
Electricity | Electricity : losses 
$ requirements 
kW kW kW % kW % kW % 
2021 | 19 762 930 | 17 682 794 | 2080135 | 10.53 | 1731668 | 8.76 | 348467 | 1.76 
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753 877 876 100 776 
2020 18 7 16 7: 2 2. 1124 1 7 8.04 ... 2.30 
2019 19 7 16 ee 2 5 11.54 1 7: 9.66 .. 1.88 
2018 18 77 16 77 2 77 13.14 1 2. 9.75 x.” 3.39 
2017 18 7: 15 7 2 2. 15.19 1 7 9.85 о 534 
2016 19 7 16 7 3 7 16.67 1 77 10.04 5 6.63 


According to experts, the amount of actual losses in 0.4 - 500КУ electric networks should 
not exceed 4-6% of the energy transmitted (received) to the network. For the Russian Federation, 
this percentage is around 7-9%. In general, the maximum actual losses should not exceed 10 - 12% 
of the electricity supplied to the network . If the actual energy losses are higher than 10 - 12%, this 
is usually explained by the presence of extreme losses. As can be seen from the tables, such a 
situation exists in many countries. [2] 


According to the voltage level of electric networks, the limiting relative technological losses 
of electric energy according to the electric energy released into the network should be as in table 7: 


Table 7 
Relative technological losses 
220 - 500 kV 2-4% 
110 kV 4 - 6% 
35 КУ 6 - 8% 
6-10КУ 8 - 10% 
0.4kV 10 - 12% 


Experience shows that most losses occur in the low voltage network, i.e. 10/6/0.4КУ 
networks and It is in those networks that the greatest difficulties in reducing losses occur. 

The reasons for this are mainly the followings: 

- lack of measuring devices or expired measuring devices of some consumers; 

- incorrect readings of balance and consumer meters; 

- intentional violation of the meter scheme and the presence of unaccounted, consumed 
(stolen) electricity; 

- methodical reporting due to the absence or failure of meters, voltage and current 
transformers, etc. 

As can be seen from the tables, in recent years, the dynamics of electricity losses in the 
electricity grids of many foreign countries have been decreasing. During the mentioned period, the 
losses in the networks of some countries decreased by about 4-5 times due to the installation of 
modern meters for electricity consumers, improvement of control of electricity consumption, data 
collection and processing, automation of the management system. 

For the Republic of Azerbaijan, these indicators are more encouraging. Thus, the losses in 
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the city of Baku have been reduced to 7-8 times in recent years, to 5%, that is, to the level of 
developed European countries. However, the actual loss rates in the regions still remain high. 


The connection between the increase of actual losses in the networks and the economic 
crisis is characteristic both for the Russian Federation and for other former Soviet Republics, which 
have entered the period of transition from centralized methods of economic management to partial 
market methods. It is clear that this is due to the weakening of control over electricity consumption 
in such a period, as well as the increase in the number of thefts of electricity, as well as the 
deepening of problems due to the decrease in the ability of a significant part of consumers, 
primarily the population, to pay. 

However, the analysis of the dynamics of energy losses both in Azerbaijan and in other 
former Soviet Republics shows that there is practically no significant reason for the increase in 
actual losses. The fact that actual losses are still high in the regions of our country can be mainly 
explained by the following reasons: 

- Slight increase in crown losses due to excessive reactive power in 110kV lines during off- 
load hours; 

- increase in the load of low-voltage networks due to the increase in the share of electricity 
consumption of household subscribers; 

- increase in non-technical (theft) losses; 

- relative increase and differentiation of wholesale and retail prices of electricity; 

- reduction of government subsidies (debt write-off) to electricity grids compared to 
previous years; 

- distribution and transmission networks are responsible for all losses of electricity; 

- transfer of low-voltage regional networks, which have not been reconstructed in the last 30 
years, to the balance of "Azerishiq" OJSC since 2015, etc. 

It is necessary to use not only the best local but also global experience to eliminate these 
causes with practical measures and reduce the actual losses in the country's energy system to 6% in 
accordance with the long-term development strategy of the energy sector for the period until 2025. 
To solve these problems, it is important to assess the actual losses of electricity in transmission and 
distribution networks, as well as to mobilize resources to reduce these losses. However, in order to 
estimate losses and reserves, it is necessary to study the stages and structure of the balance of actual 
losses of electricity, including non-technical (commercial) losses and the factors affecting them in 
as much detail as possible. It is proposed to establish the stages of electricity balance in distribution 
networks according to Figure 1: 


Figure 1. Stages of electricity balance in distribution networks 
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The balance of electricity by voltage steps should be drawn up in distribution networks as follow 


The balance of high-voltage power networks is 110-35 kV networks 


Balance on low-voltage electric networks 10/6/0.4 kV networks 


Energy system balance sheets of countries: 


Overall energy system balance. Production-transmission-distribution-sale 
Balance of 500/330/220/110 kV networks by transmission system operator 


Balance on distribution networks 110/35/10/6/0.4 kV networks 


The state of the country's energy system is mainly characterized by a single fuel-energy 
balance model that combines the production and consumption balances of individual energy 
carriers. According to the definition [2]: "The fuel energy balance is the ratio of the amount of fuel- 
energy resources received as a result of production or import for an economic object or a certain 
area, and the decrease due to local consumption or export." Thus, the main task of the fuel-energy 
balance is to show the real structure of production, transmission, distribution and use of energy 
resources in a certain area. The fuel-energy balance can be compiled at the state level or at the level 
of its separate territory (even enterprises, settlements). 


Another important problem is the prediction of the fuel-energy balance. In the Republic of 
Azerbaijan, the fuel-energy balance is prepared and approved by the Cabinet of Ministers for each 
year. This document also provides forecasts for the next 3 years. We believe that there is a need to 
further improve the design and forecast scheme of some components of the fuel-energy balance and 
to further strengthen the control over its implementation. As a result of unclear or non-transparent 
fuel-energy balance, the following problems may arise in energy policy [2]: 

- reduction of technical and economic availability of energy services; 

- low economic efficiency due to high losses in the fuel-energy complex (mainly distribution 
networks); 

- decrease in energy security of the country; 

- slowdown of economic growth; 

- acceleration of inflation; 

- the increase in the payment burden in the family budget of energy carriers and the decrease 
of payments; 

- increase in the burden of utility payments and increase in receivables in city and district 
budgets. 


In order to develop a high-quality interconnected fuel energy balance and forecasts, 
systematically integrated balances should be drawn up at all levels, from bottom to top, in all areas 
and enterprises that are components of the balance. To solve this global problem, divide the task 
into parts is appropriate. 
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I believe that in the preparation of high-quality electricity balance in electric networks, 
attention should be paid to 3 problems: 
1. Absence of reliable and centrally controlled full accounting of electricity received 
(received), outgoing (transmitted to other networks) and consumed (certified software); 
2. Absence of the legislative framework for the formation of the balance in electric 
networks, as well as the absence of approved legal normative documents-rules; 
3. The need to improve the normative-methodological base and terminology used in the 
formation of the structural components of the electric energy balance and tariff regulation. 
In addition, terminology and definitions related to losses and sales of electricity used in 
management need to be clarified. Thus, there are still ambiguous interpretations in various 
documents, separate scientific articles, and documents of economic subjects. The lack of an 
approved terminology in this field prevents integration between separate structures of the economy. 
Professionals find it difficult to fill out forms and make mistakes. 


I should mention that back in 2004, in OJSC "Bakielektrikshabek" (at that time it was 
managed by the Turkish company "Barmek Holding"), for the first time, a proposal was made by 
me to calculate the balance of the electricity consumed in the building by installing a balance meter 
in the distribution box at the entrance of multi-storey buildings. This was aimed at concretely 
finding the location of electricity thefts in apartments and acting accordingly. Although there were 
skeptical approaches to this proposal by Turkish experts, with the help of our local experts, balance 
reports were drawn up for multi-storey buildings, and later for TM and KTMs, software was 
developed, and for the first time in our republic, a team of 3 people at "Bakielektrikshabek" OJSC 
made a balance and a loss investigation department was established. Today, about 200 specialists 
work in the Purchase and Balance Department at "Azerishiq" OJSC, balance reports are drawn up 
for regions, cities, districts, areas, TM and KTMs, high and low voltage networks. The 3 problems I 
mentioned above create difficulties in the more perfect preparation of balance sheets. 


I should also mention that the balance sheets of the networks may not coincide with the 
division of the administrative-territorial units of those cities and regions. The main issue here is the 
historically formed form-design of networks. The main criterion is that the network balance, at any 
stage, should be complete, allowing the calculation of network losses, their localization and 
determination of their nature (type), and the electricity received and transmitted to other networks 
can be accurately measured . Therefore, it is possible to attribute some villages belonging to another 
region to the network of the other region or vice versa. Taking into account the complexity of the 
networks of some districts and cities, it is more efficient to divide them into several networks and to 
conduct separate balance sheets. 


For example, Bilajari ESH, which is part of Binagadi district of Baku city, has been operating 
as a separate network since 2004 and a separate electricity balance is drawn up. Also, Zira ESH of 
Khazar region and Mashtaga ESH of Sabunchu region operate independently. Baglar ES, which has 
no administrative-territorial unit, has been operating separately since 2006 and is almost the most 
difficult and loss-making farm of the electricity networks around Baku. Designing the networks in 
this way, drawing up the balance sheets had a great impact on the reduction of the losses of the 
electricity industry of those regions, the Absheron Peninsula in general. 
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The remainder of the actual loss 
(portion exceeding the norm) 


. Normative Commercial Other non- 
Technical losses losses of losses technical 


Load losses Losses due to weather Conditionally stable 


Figure 2. The structure of actual losses of electricity in electrical networks. 


Technical losses 


Load losses Conditionally stable losses Climatic losses 
In power transformers In power transformers Crown losses 


In overhead and cable lines In transformers of arc choke In leakage currents in 
reactors insulators of overhead lines 


In support and other 
elements of overhead lines 


Normative technological losses (waste) 


£ 
wa 


In compensating devices 
(SK, SKB, STK) and 
compensating devices with 


Loss of ice melting 


In high frequency 


ecammiuinieatian davirac 


In surge arresters and 
switch-type arresters 


b) In current choke reactors 


In Y/S and connections of 


busbars and wires in PQ In the connectors of high- 
frequency communication 
In current transformers - : ; 
In the insulation of cable lines 
In single-ended meters 


In voltage transformers 


Other losses due to 
negligence 


All the indicated practical measures made it possible to shape the balance sheet of 
"Bakielektrikshabek" OJSC starting from 2004, and "Azerishiq" OJSC since 2015, to minimize 
losses. However, the lack of an official Methodical instruction registered in the Ministry of Justice 
for the formation of structural components of electricity balances in the country does not allow for a 
fully reliable assessment of the scale of electricity losses and reduction of electricity losses. This, in 
turn, complicates the development of feasibility studies and programs for reducing losses during the 
period of tariff adjustment and in the long term. 
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Let's consider the structure of the actual losses of electricity in electrical networks 


Technological losses (waste) 


Technological losses during distribution 


Non-technical (commercial) losses 


Technological losses during the sale 


Errors of measuring 


с) 


Losses due to weather 


Errors of measuring devices 


Thefts 


Failure to simultaneously remove indicators of measuring 


Missing or incorrect links system 


Subscribers whose meter readings are not read 


Subscribers whose meter readings are not fully 


Unregistered users 


Software errors 


Debts (past due) 


Losses due to negligence 


Figure 3. Structure of technical (a), technological (b) and commercial (c) losses 


The structure of technical losses of networks has been studied almost enough in our country 
and substantial scientific researches have been conducted. Figure 3, a, b, c shows the structure of 
technical, technological and commercial races: 


The above tables (schemes) are a traditional approach to actual losses. In simple terms, the 
actual losses in the distribution system are divided into technical and non-technical. Technical 
losses are related to energy consumed in conductors and equipment used for energy distribution, 
depending on the technical characteristics of the distribution network. These losses can occur 
anywhere in the electrical distribution system. Technical losses in power distribution systems can be 
divided into losses depending on the intensity of the current in the elements of the electric circuit or 
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losses independent of the intensity of the current. According to the load condition, technical losses 
are divided into two parts, depending on the load and independent of the load [9]. Technical loss 
due to electromagnetic induction also occurs in the core of transformers. There are many methods in 
the modern literature for accurate calculation of technical losses, even special software is operating 
at the Azerbaijan Scientific Research and Project Research Energy Institute. [5] 


For example, Bilajari ESH, which is part of Binagadi district of Baku city, has been operating 
as a separate network since 2004 and a separate electricity balance is drawn up. Also, Zira ESH of 
Khazar region and Mashtaga ESH of Sabunchu region operate independently. Bağlar EŞ, which has 
no administrative-territorial unit, has been operating separately since 2006 and is almost the most 
difficult and loss-making farm of the electricity networks around Baku. Designing the networks in 
this way, drawing up the balance sheets had a great impact on the reduction of the losses of the 
electricity industry of those regions, the Absheron Peninsula in general. 


All the indicated practical measures made it possible to shape the balance sheet of 
"Bakielektrikshabek" OJSC starting from 2004, and "Azerishiq" OJSC since 2015, to minimize 
losses. However, the lack of an official Methodical instruction registered in the Ministry of Justice 
for the formation of structural components of electricity balances in the country does not allow for a 
fully reliable assessment of the scale of electricity losses and reduction of electricity losses. This, in 
turn, complicates the development of feasibility studies and programs for reducing losses during the 
period of tariff adjustment and in the long term. 


Let's consider the structure of the actual losses of electricity in electrical networks 


Causes of non-technical losses 


Immeasurable energy Payment errors Payment losses itkilər 
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Illegal connections Unpaid bills 
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Faulty connections Registration of adjusted Payment losses 
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Inconsistency in measuring Software errors 


devices 
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Consumers without a code of follow-up 


Malfunctions and errors in the 


Errors in reading indicators 
software 


Figure 4. Information on the causes of non-technical losses 
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Non-technical losses can be defined as the difference between the consumed electricity and 
the measured electricity or billing amount. We call it commercial losses. Studies of these losses 
show that the ratio of non-technical losses to total distribution losses cannot be accurately 
estimated. Commercial losses depend on the technical condition of the networks, the social welfare 
of the urban population of the region, the organization of network management, the compliance of 
the staff serving the population, etc. it varies depending on other factors. However, commercial 
losses can be estimated to be 20-30% of the actual losses, and in some cases up to 50%. 


A separate article will be presented on the structure, causes and elimination methods of non- 
technical losses. Currently, we have tried to summarize the information about the causes of these 
losses in the following table (Figure 4) [8]. 


Electricity is an indispensable element of our daily life. As a result of incorrect and 
inefficient use of this energy, technical and non-technical losses occur at the stages from 
production, that is, from power plants to the point of consumption. Reducing these losses can allow 
more efficient use of existing energy systems and thus maintain a balance of demand and supply in 
the electricity sector. 


CONCLUSION. 


1. Energy losses in transmission and distribution networks around the world and their 
monitoring and gradual reduction are becoming more important issues. In a sense, efforts to reduce 
energy losses are considered equivalent to research to produce energy. 

2. The investigation of reporting materials and statistical data from foreign sources shows that 
the main criterion of the economic efficiency of the production and distribution of electricity is the 
ratio of the actual losses of electricity to the electricity produced by the generators of power plants. 

3. Technical losses, which are a part of actual losses, have been sufficiently researched in our 
country, relevant scientific works have been written, methods of calculating technical losses have 
been investigated, and regulations for electric networks have been developed. However, it is 
impossible to accurately calculate non-technical losses by mathematical methods. Application of 
balance sheets in transmission and distribution networks, specification of its structure and 
components, allows determination of actual losses, as a result, non-technical losses. 

4. Carrying out balance reports from A to Z on all stages allows to localize the actual losses, 
predict the presence of commercial losses in specific 10/6/0.4kV lines and MVs and take measures 
to reduce them. 

5. There is a need to study non-technical losses in distribution networks in our country, to 
investigate the causes of their occurrence, ways to reduce them, to study the effectiveness of 
methods of preventing electricity theft, to study the relationship between the actual losses and the 
dynamics, and to make concrete proposals. 
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ABSTRACT. 


The article discusses the features of a mathematical model developed to study the processes 
in the system of a two-cylinder hermetic piston compressor when it operates in an open cycle and 
refrigerant - 22. Some results of computational experiments are given that are in good agreement 
with the literature data. 

Keywords: Hermetic piston compressor; mathematical model; process modeling; calcula- 
tion of flow parameters; explicit difference scheme. 


INTRODUCTION. 


As is known, in order to create refrigeration equipment that meets the requirements of world 
standards, it is necessary to increase the level of research work using the latest equipment, research 
methods and processing of their results. 


At the Azerbaijan Technical University, a computer method for conducting a mathematical 
experiment on a set of processes occurring in a refrigeration unit of a household air conditioner, 
which has a rotary compressor and operates in cold and heat modes, has been developed and 
implemented; in dynamics, including an improved mathematical model that takes into account the 
unsteady wave nature of the refrigerant flow in pipelines, hydraulic resistance, variability of 
physical characteristics of the medium, phase transformations, heat and mass transfer, the influence 
of moving elements, narrowing and expansion of channels. The model allows for computational 
research in time and coordinates, and thereby turns the system, as it were, into a “transparent” one 


[1]. 


Due to the fact that compressors are widely used not only in air conditioning systems, but also 
in other refrigeration systems, and they can also be used as compressor units for various purposes, 
we have developed a methodology that allows us to conduct computational studies of flow 
parameters in small sealed volumetric systems compressors when operating in an open cycle. In this 
case, the model simulates the operation of a compressor on a stand to determine air capacity [2]. 
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Studies of the operation of small hermetic compressors (XqRV - 1.75 and XKV - 6) have shown 
that this technique is quite effective. 


The research and accumulated experience made it possible to develop and implement on a 
modern computer a method for calculating the parameters of processes occurring in the system of a 
hermetic two-cylinder piston compressor PG - 5 and operating on freon - 22. Technical data and 
characteristics of the compressor are given in [4]. 


The choice of such a compressor is due to the fact that it has a relatively large cooling 
capacity (5.8 kW in “standard” mode), the flow part is similar to the flow part of medium and large 
compressors, and the multi-cylinder nature allows one to take into account the displacement of 
processes depending on the angle between the shaft elbows. 


The design diagram of the compressor is shown in Fig.1. 
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Figure 1. Compressor design diagram. 
1 - cavity of the sealed casing; 2.5 — supravalvular cavities of the suction 
and discharge valves; 3, 4 — cylinder cavities (suction and compression); 
6 — receiver; I, II, Ш- pipelines. 


In the developed model, the pressure in the pipelines of the system is determined from the 
general equations of gas dynamics (equations of motion and continuity) presented in model form 


[3]. 


04 84 

ЕЕ vo fw) (1) 
06 ol 

2.50 (2) 


VVhere, q=@O+w: {= Ф - у; b=a+w (3) 
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x. 
c=a-w; f(w) as Pe 


pcd (4) 
ap 


The temperatures of the refrigerant flow and the walls of the pipelines were determined from 
the energy and heat balance equations 


öl бі 
pwf — + of a = ash; (Ber -T) (5) 
öx ör 
д0ст 
Ayhy (Т, 022) 8 n Cm öz = ah, (Ber -T) (6) 


Qi jii 7 Wk j di; ^ dii; 
аа сеа £(w,;) (7) 
Pa — 27 aL 
ee ue UN (9 
Р-Р. 
Pa ии (9) 
к.) 


і і, AES E 
( wu oath, = (eta) On Taa) (10) 


0, 


ie. 
8 "Ст = = (aha), ‚©, Tys) (eh; hh. -в, (11) 


The stability of calculations and the convergence of solutions were ensured by compliance 
with the Courant condition: 


(12) 


In relations (1) = (12), W, m/s - is the flow speed; , Ф, m/s - substitute function for bringing 
the gas dynamics equations into the form of model equations; , 4, m/s - speed of sound; b, c - the 
speed of propagation of the pressure pulse through the pipeline in the forward and reverse 
directions, respectively; P, Pa - pressure of the working fluid in the flow; , p , kg/m3 - density of the 
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working fluid; R, m - internal radius of the pipe; б . coefficient of hydraulic resistance; Í , C/kg - 
enthalpy; 75 Wt/m2 К - heat transfer coefficient from the channel wall to the refrigerant flow; 
O5. Wtm2 К - heat transfer coefficient from the blown medium to the channel wall; , К - 


hy ,h; 


temperatures of the channel wall and working fluid; , m2/m - specific outer and inner 


surfaces of the channel; 8% , kg/m - mass of channel material per unit channel length; бт, C/kg К - 


specific mass heat capacity of the channel material; r(sek) and x (m) - coordinate and time; and - 
numbers of steps along the pipeline length and in time, respectively. 


To determine the parameters in the cavities of the system, equations of boundary conditions 
for the compressor and for various cavities were used. 


The equations for the boundary conditions of the compressor consist of the equation of 


motion of the piston, changes in the volumes of the suction and compression cavities, the equations 
of motion of the suction and discharge valves, etc. These equations are composed similarly to the 


equations given in the work Ш. 

When determining the parameters of the working fluid in the cavities of the system, the 
equation of the first law of thermodynamics for a body of variable mass, the mass balance equation 
and the equation of state were also used. 

The equation of the first law of thermodynamics in differential form is: 


dQ + Е, = dU +dL+dE, (13) 


Mass balance equation 
dM = dM,, — ами, (14) 


Equation of state (Bogolyubov-Mayer equation): 


r 5р. | 
а-ы 2; as 


т Tj 


Т 
Where, Z= p-107/RƏT, Е = 8,31437/ u; т--- (16) 


kp 


In expressions (13) “ (16) dQ. heat imparted to the working fluid from the outside; dE, _ 


the amount of thermal energy entering the working space (through the outlet valve, leakage through 
cracks), 40 - change in the internal energy of the working fluid; dL . external work on the working 
fluid; dE, the amount of energy lost from the working space (through the discharge valve and 


leakage through the cracks), dM . change in the mass of the refrigerant in the working cavity; 
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GM, . mass of the working fluid entering the working cavity; dM... mass of the working fluid 


displaced from the working cavity; T, K - temperature; P, Pa - pressure; M. molecular mass; P, 


g/sm3 - density; Tkr, К - critical temperature; “/- virial coefficient depending on temperature. 


The individual components of equations (13) and (14) are calculated using similar equations 
given in work [1], the virial coefficients of equation. (15), as well as thermodynamic, 
thermophysical and heat - mass transfer characteristics and gasdynamic resistances are determined 
using the formulas given in [5]. 


The final calculated parameter in the model is the volume flow through the nozzle in the 
receiver (Fig. 1). The diameter of the nozzle is determined from the condition of equality of the 
actual volumetric cooling capacity of the compressor and the flow rate through the nozzle. 

The mass flow rate through the nozzle per cycle is determined according to 


G- Уи №. р,Ат, kq/sec (17) 


i=l 
and the actual volumetric capacity of the compressor 


RT 
V. 7.7 m?/sec (18) 


t. P; 


Where, 77 f. effective flow area of the nozzle, m2; 
We - flow velocity leaving the receiver, m/s; 


Pr. density of the working fluid in the receiver, kg/m3; 


Polo. pressure and temperature of the medium into vvhich the 


refrigerant flovv exits, Pa, K, 
R - gas constant of the vvorking fluid, C/kg K, 
AT - time step, sec; 
7,- time of one cycle, sec. 


The program was compiled in FORTRAN-e, the calculations were made on a personal 
computer. 


Figure 2 shows the calculated indicator diagram of the compressor when it is operating in 


—. 0 — 0 
the “standard” mode (Io 5. 


in the sealed casing, in the over-valve cavities of the suction and discharge valves and in the 


). The diagram also shows graphs of pressure changes 


receiver [6]. 
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Figure 2. Compressor design indicator diagram. 
Rvp, Rvk, Rnk, Pres - pressure in a sealed casing, in the over-valve 
cavities suction and discharge valves and in the receiver; 1 - traffic 
schedule suction and discharge valve plates. 


As you can see, the indicator chart has a theoretically imaginable form. However, the 
opening of the suction and discharge valves is delayed, and the closing of the discharge valve 
occurs during the suction process of the next cycle. Apparently, this is due to the peculiarity of the 
crank mechanism (uneven movement of the piston) and the low values of the polytropic expansion 
and compression index for freon - 22, which practically brings these processes closer to isothermal. 
It should be noted that by optimizing the flow sections of valves and spring parameters, it is not 
difficult to shift the opening and closing moments of valves in the desired direction. 

Figure 3 shows a graph of pressure changes along the III tube (Fig.1), connecting the 
compressor cavities to the receiver. As can be seen, the intensity of the pressure wave is not great. 
We carried out calculations of the actual volumetric productivity of the compressor in various 


ty — -30--0”C 


modes (vvithin and 1 = 30 $55'€3 and determined the values of the compressor 


supply coefficient [7]. 
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Figure 3. The change in pressure along the III tube at 
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Figure 4 shows graphs of changes in the supply coefficient depending on the boiling and 
condensation temperature. The results obtained are in good agreement with the literature data [3]. 
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Figure 4. Changing the feed ratio depending on the modes compressor operation. 
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XULASO 


İşdə (p.p. T) asılılıqlarının ölçülməsi zamanı hər bir müvazinət halında atmosfer təzyiqində 
qrafik ekstrapolyasiya yolu ilə sıxlığın yüksək dəqiqlikli qiymətlərinin alınması məqsədi ilə 
təzyiqin maksimum dərəcədə mümkün olan aşağı qiymətləri yaradılmağa səy göstərilmiş, alınmış 
qiymətlər sıxlığın DMA 5000M qurgusunda ölçülmüş qiymətləri ilə müqayisə olunmuşdur. 
Azərbaycannın Xaçmaz rayonunun “4-cü bölmə” termal suyunun müxtəlif təzyiq və temperatur- 
larda sıxlığının təcrübi qiymətləri ölçülmüşdür. Öncə təcrübə qurğusunun ölçmə dəqiqliyinin 
yoxlanılması üçün su, toluol üçün alınmış nəticələrin müxtəlif ədəbiyyatlarda verilmiş məlumat- 
larla müqayisəsi aparılmışdır. Alınmış təcrübi nəticələr cədvəldə göstərilmiş və hal tənliyi ilə 
ifadə olunmuşdur. 


Açar sözlər: sıxlıq, təzyiq, temperatur,termal sular. 
SUMMARY 


During the measurement of (p, p, T) dependences, in order to obtain high-accuracy values of 
density by means of graphical extrapolation at atmospheric pressure in each case of equilibrium, 
efforts vvere made to create the lovvest possible values of pressure, the obtained values vvere 
compared vvith the values of density measured on the DMA 5000M device. compared. The 
experimental values of the density of the "4th division" thermal vvater of Khachmaz region of 
Azerbaijan at different pressures and temperatures were measured. First, in order to check the 
measurement accuracy of the experimental device, the results obtained for water and toluene were 
compared with the information given in various literature. The obtained experimental results are 
shown in the table and expressed by the equation of state. 


Keywords: density, pressure, temperature, thermal waters. 
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Currently, extensive measures are being taken both within the Republic and internationally 
to reduce the amount of carbon dioxide emitted into the atmosphere. Great work is being done in 
this direction in Azerbaijan (including at the Department of Energy Efficiency and Green Energy 
Technologies of Azerbaijan Technical University). In accordance with the "State Program on the 
Use of Alternative and Renewable Energy Sources in the Republic of Azerbaijan" approved by the 


Decree of the President of the Republic of Azerbaijan dated October 21, 2004, the Ministry 
of Industry and Energy of the Program the implementation coordinator has been determined [1]. In 
the State Program, the main directions for the implementation of measures related to the use of the 
potential of wind energy, as well as solar energy, the hydropower of geothermal waters, mountain 
rivers and water channels, as well as the energy of biomass, as the most efficient energy sources for 
our Republic, have been defined. 


In addition to the regional fracture in the zone of thermal water distribution, there is also a 
zone of tectonically disturbed rocks with a branched system of sharply falling (descending) cracks 
covering the entire Mesozoic layer. The velocity of water is higher in the main drainage channels, 
i.e., where the transfer of water heat to the surrounding rocks is the least, thus the water has the 
maximum temperature, and in the zones with many cracks. This is confirmed by the outputs of 
geothermal sources in the erosion depressions of the intersection of the main sututar zones [2]. 


The hydrogeological parameters of the discharge zone of thermal waters directly depend on 
the degree of cracking and the character of the hydrous rocks. The abundance and abundance of 
gryphons in natural outcrops indicates a well-developed rift. 


Geothermal energy resources of the studied region (Khachmaz region) are particularly 
valuable. Thus, the colorful chemical composition, high healing properties, favorable geographical 
position of the area create a good basis for their treatment purposes and wide application in various 
fields of the national economy [2,3]. 


The analysis of the chemical composition of thermal water of Khachmaz district "4th 
section" was measured in atomic emission spectrometer with IRIS Intrepid II Optical Emission 
Cheomotograph inductively coupled plasma [3,4]. The results show that the majority of chemical 
elements are sodium (Na). Na constitutes approximately 72.41 + 90.12% of all chemical substances 
in the thermal water of the Khachmaz region of Azerbaijan. In the tables 1-4 shown below, the 
geographical coordinates of the thermal water of the "4th section" station located in Khachmaz 
region of Azerbaijan, the temperature at the time of exit from the source and the amount of minerals 
in the chemical composition are given [3,4] 


Geographical coordinates of Khachmaz district "4th section" thermal water and its 
temperature at the moment of exit from the source 


Table 1. 
Name of source Geographical coordinates Temperature at exit 
North 41736/17” 
“Unit 4” 5: Т = 327,15 К 
East 48 41'54" 


Chemical composition of thermal water of "4th division" of Khachmaz region of Azerbaijan 
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Table 2. 
The А11670 | As1890 | B2089 | Ba2304 | Ca3181 | Cd22 | Co2286 E1205 İCu32 Fe25 
amount 88 5 47 199 
of 
minerals 
in the 
sample 
Section | «0.01 «0.01 1.87 0.11 158.0 | «0.001 | «0.001 | <0.0 | «0.0 | 0.0 
4, 1 1 2 
mg/liter 

Table 3. 
The Hg184 | K766 | Li670 | Mg279 | Mn293 | Мо204 | Na81 | Ni231 | P2136 
amount of 9 4 7 0 9 5 83 6 
minerals in 
the sample 
Section 4, | <0.02 | 18.30 | 0.16 28.20 0.11 0.02 832 | <0.01 | <0.01 
mg/liter 


Table 4. 
The Pb2203 | $1820 | Sb217  Se196 | $12124 | Sr407 | Ti334 | T1190 | V2924 
amount of 5 0 7 9 8 
minerals in 
the sample 


Section 4, <0.01 | 37.40 | <0.02 | <0.02 3.09 9.28 | <0.01 | <0,05 <0.01 
mg/liter 


Taking into account the fact that the vibrating tube densimeter device needs to be calibrated 
with at least two substances for the study of the density of liquids after checking the device's 
performance, for this purpose, water, toluene are used as the main calibrator. standard) were 
selected as substances. The results obtained for water and toluene were compared with the 
information given in the literature. As a result of the comparison, the difference between the 
obtained values for the density of water and toluene and the information in the literature shows that 
the estimated errors of the measurements in the device are very small. Obtaining results with a 
small error and close to each other shows the high accuracy of the created experimental device 
[4,5]. 


After the verification experiments on the properties of water and toluene (p,p,T) were 
carried out, the density of the thermal water of Khachmaz district "4th section" of Azerbaijan was 
measured at high pressure and different temperatures in the experimental facility that works with 
the vibrating tube densimeter method. . During the measurement of (p,p,T) dependences, in each 
case of equilibrium, an effort was made to create the lowest possible values of the pressure in order 
to obtain high-precision values of the density by means of graphic extrapolation in the atmospheric 
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pressure, and the obtained values are the DMA of the density It was compared with the values 
measured in the 5000M device. The values obtained by different methods agree well within + 
0.02% [3-5]. Researches for thermal water of Khachmaz district "4th division" were conducted at 
temperatures T=(278.15+373.15) K and pressures up to p=40 MPa. Experimental indicators 
obtained on (p,p,T) dependencies are given in table 5. 


Experimental values of the density of the "4th division" thermal water of Khachmaz region of 
Azerbaijan at different pressures and temperatures 


Table 5. 
P ЕЕ > T P. Pe 3 
MPa kql m? K MPa kqlm” K 
0.201 1007.44 278.15 0.214 991.53 328.15 
5.006 1009.65 278.15 5.006 993.57 328.17 
10.006 1011.92 278.16 10.301 995.81 328.16 
15.214 1014.26 278.14 15.921 998.17 328.15 
20.004 1016.39 278.15 20.152 999.93 328.14 
25.301 1018.72 278.16 25.008 1001.94 328.15 
29.986 1020.76 278.15 30.102 1004.03 328.13 
35.114 1022.97 278.13 35.026 1006.03 328.15 
39.997 1025.04 278.15 39.996 1008.03 328.15 
0.690 1006.06 288.16 0.302 983.52 343.15 
5.061 1007.95 288.11 5.014 985.64 343.17 
10.162 1009.96 288.10 10.006 987.85 343.15 
15.166 1011.85 288.13 15.308 990.14 343.16 
19.942 1014.07 288.11 20.410 992.30 343.15 
24.996 1016.13 288.13 25.008 994.21 343.16 
30.010 1018.40 288.11 29.998 996.23 343.16 
34.744 1020.45 288.13 35.047 982.30 343.15 
40.002 1022.67 288.13 39.995 1000.14 343.15 
0.798 1002.95 298.19 0.304 974.18 358.15 
5.132 1004.77 298.19 5.008 976.40 358.14 
9.979 1006.82 298.19 10.009 978.69 358.16 


14.879 1009.03 298.14 15.308 981.04 358.15 
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20.098 
25.123 
30.022 
34.932 
39.846 
0.385 

5.158 

10.079 
15.112 
19.962 
25.287 
29,992 
35.037 
40.130 


1011.23 
1013.33 
1015.38 
1017.28 
1019.24 
997.61 

999.87 

1002.12 
1004.32 
1006.44 
1008.74 
1010.62 
1012.69 
1014.69 


298.13 20.008 
298.13 25.021 
298.12 29.987 
298.13 35.030 
298.14 39.994 
313.07 1.393 

313.08 5.416 

313.08 10.407 
313.08 15.587 
313.06 0.420 

313.06 25.274 
313.07 30.049 
313.06 35.201 
313.08 39.745 


983.07 358.15 
985.18 358.16 
987.20 358.17 
989.19 358.15 
991.08 335.15 
963.977 373.08 
965.99 373.08 
968.36 373.08 
970.77 373.09 
972.95 373.09 
974.97 373.09 
977.01 373.09 
978.99 373.09 
980.72 373.10 


The obtained experimental results are expressed by the following equation of state: 


р=Ар'+Вр'+Ср” (1) 
Here, coefficients A(T), B(T) and C(T) depend on temperature in polynomial form: 
3 2 , 2 | 
A(T)= Ха", ВОТ) = Xe", CT) = Хет (2) 
Values of coefficients а, b and с in equation (2) are given in table 6. 
Tablel 6. 


a1 = -1.2533813 


bo = 736.75489731 


Co= -488.2859137 


a = -0.339952- 107 


bı = -3.43166265 


сі- 3.2130145 


аз = 0.96906297- 107 


bə = 0.87206243- 107 


c2 0.61334197-10? 


The equation (1) allovvs to express the experimental values of the dependence of thermal vvater 
(р,р,Т) of "Section 4" with average error of 0.007%, taking into account the values of the 
coefficients A(T), B(T) and C(T). gives 
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